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Simulating  the  Management  of  Even-Aged  Timber  Stands 


Clifford  A,  Myers 


Introduction 

Good  decisions  that  lead  to  efficient,  low- 
cost  operations  are  as  necessary  for  timber 
management  as  for  any  other  business  activity. 
Gathering  of  information  on  which  to  base  such 
decisions  is  unusually  difficult  because  of  the 
peculiar  characteristics  of  the  business  of  grow- 
ing wood.  The  time  needed  to  mature  a  crop 
is  long,  and,  with  many  tree  species,  results  of 
changes  appear  slowly.  Decisions  must  be  made 
for  large  areas  of  land  and  for  numerous  pos- 
sible combinations  of  biological  and  economic 
conditions.  A  procedure  whereby  a  manager 
can  examine  quickly  the  probably  outcomes  of 
many  variations  in  management  should  be  of 
value.  Such  a  procedure  appears  especially 
useful  when  the  computation  time  required  is 
about  10  seconds  for  each  combination  of  pos- 
sible alternatives  examined.  This  approach  to 
managerial  decisionmaking  is  made  possible 
by  electronic  computers  and  the  technique  of 
simulation. 

Simulation  involves  the  creation  and  opera- 
tion of  a  model  that  is  similar  in  relationships 
to  the  natural  system  studied.  The  system  may 
be  a  factory,  a  forest,  or  any  one  of  countless 
other  subjects  of  interest.  Mathematical  simu- 
lation, where  the  model  is  a  set  of  equations, 
is  a  powerful  tool  of  management  and  investiga- 
tion. Mathematical  models  are  used  to  discover 
new  facts  and  to  test  alternatives.  Time  can 
be  compressed  so  a  sequence  of  events  can  be 
studied  in  a  fraction  of  the  time  required  for 
operation  of  the  natural  system,  thereby  permit- 
ting projections  into  the  future (Choraf as  1965). 

Manipulation  of  a  model  by  simulation  is 
especially  useful  where:  (1)  the  actual  system 
should  not  be  changed  until  decisions  have 
been  made;  (2)  possible  changes  are  numerous 
and  are  costly  in  money,  effort,  and  time;  and 
(3)  exact  data  cannot  be  obtained  because  of 
cost  or  complexity  of  the  relationships  involved 
(Churchman  et  al.  1957,  Dane  1966). 


Several  authors  have  presented  the  prin- 
ciples of  simulation  to  foresters,  and  have  shown 
how  forestry  simulators  may  be  constructed. 
These  authors  have  discussed  the  development 
of  models  (Chappelle  1966)  and  the  systems 
approach  to  forestry  problems  (O'Regan  et  al. 
1966).  Others  have  described  simulators  of 
timber  production  activities,  and  presented  in- 
formation obtained  from  their  operation  (Clut- 
ter and  Bamping  1966,  Gould  and  O'Regan 
1965).  Thompson's  (1966)  discussion  of  the 
need  for  tests  of  alternatives  in  forest  regulation 
points  out  a  value  of  simulation  for  research. 
Validity  of  many  long-accepted  concepts  of 
timber  management  can  be  tested  quickly, 
cheaply,  and  for  a  wide  variety  of  conditions. 

Gould  and  O'Regan  (1965)  performed  a 
special  service  to  timber  managers  by  publishing 
a  simulator  at  a  time  when  examples  were 
relatively  inaccessible  to  most  foresters.  Their 
computer  program  will  long  serve  as  a  source 
of  ideas  and  procedures;  it  has  so  served  for 
the  program  presented  here. 

Program  MANAGD  (appendix  1)  is  a  means 
for  simulation  of  the  management  of  even-aged 
stands  for  roundwood  and  saw  logs.  It  contains 
provisions  for  stand  growth,  thinnings,  harvest 
cuts,  planting  of  nonstocked  areas,  and  other 
changes  in  forest  conditions.  Inputs  to  the 
program  permit  wide  choice  in  the  management 
alternatives  and  stand  conditions  to  be  ex- 
amined. Possible  options  and  alternatives  are 
presented  in  the  appropriate  sections  of  the 
program  description. 

The  program  was  written  in  A.S.A.  Fortran 
IV  (X3.9-1966)  and  tested  on  a  CDC  6400  compu- 
ter. 2  Program  organization  will  permit  ready 
modification  to  fit  local  tree  species  and  com- 
puting equipment.  Once  the  program  is  in  use, 
improved  or  additional  data  and  procedures  can 
be  inserted  easily,  as  they  become  available. 

1  Trade  names  and  company  names  are  used 
for  the  benefit  of  the  reader 3and  do  not  imply 
endorsement  or  preferential  treatment  by  the 
U.  S.  Department  of  Agriculture. 
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Description  of  Program  MANAGD 

Program  MANAGD  consists  of  a  main  pro- 
gram and  11  subroutine  subprograms  (appendix 
1).  Content  and  purpose  of  each  routine  are 
given  in  the  sections  that  follow.  Definitions 
of  variable  names  are  listed  with  the  source 
program  in  appendix  1  and  in  the  listing  of 
contents  of  the  data  deck.  The  test  problem 
described  on  page  14  and  reported  in  appendix 
2  provides  additional  explanation  of  the  program. 

The  terms  batch,  test,  and  game  are  used  to 
identify  individual  simulations  performed  with 
various  groupings  of  alternatives.  The  BATCH 
name  identifies  the  entire  group  of  tests  and 
games  to  be  completed  as  a  single  job  by  a 
computer.  A  test  consists  of  one  or  more 
games,  all  of  which  are  based  on  a  single  yield 
table  and  one  set  of  stumpage  prices.  Games 
of  a  test  differ  from  one  another  in  distribution 
of  acres  by  age  classes,  area  planted,  and  limi- 
tations on  the  annual  cut. 


Main  Program 

The  main  program  calls  10  subroutines  in 
proper  sequence,  and  uses  counter  UK  to  call 
an  eleventh  subroutine  (OUTPUT1)  at  specified 
intervals.  The  first  four  subroutines  called 
read  the  data  deck,  compute  and  print  a  yield 
table,  and  compute  and  print  potential  volumes 
per  acre  at  each  year  of  stand  age.  A  fifth 
subroutine  generates  a  working  circle  with  the 
specified  number  of  acres  in  each  age  class. 
Subroutines  YEARS  and  ANUAL  are  the 
dynamic  components  of  the  program.  They 
provide  for  stand  growth,  thinnings,  harvests, 
losses,  and  other  changes  in  volume  and  value. 
Four  subroutines  produce  a  record  of  results. 
Numbers  of  acres  in  each  age  class  are  printed 
at  the  end  of  the  first  year  of  a  game  and  at 
the  end  of  each  decade.  Volume  and  dollar 
values  are  printed  at  the  end  of  each  game. 
An  optional  subroutine  (SUMRY)  prints  selected 
values  from  several  games  of  a  test  together 
to  simplify  comparisons. 

The  main  program  enters  BATCH  name  and 
the  number  of  tests  in  a  batch  into  computer 
memory. 


Subroutine  INPUT1 

INPUT1  is  called  once  each  test  to  read 
values  that  apply  to  all  games  of  a  test.  Values 


entered  include  stumpage  prices,  minimum 
commercial  volumes,  and  items  used  to  compute 
a  yield  table.  Controls  on  the  program  are 
entered  as  number  of  games  in  a  test,  number 
of  years  in  a  game,  and  number  of  columns 
of  OUTPUT2  to  be  summarized  by  SUMRY. 

Values  read  by  subroutine  INPUT1  are 
printed  as  part  of  the  output  of  other  sub- 
routines. 


Subroutine  YIELD 

Subroutine  YIELD  reads  four  tables  and 
computes  and  prints  yield  tables  for  managed, 
even-aged  stands.  It  is  called  once  each  test 
to  produce  the  yield  table  that  will  apply  to 
all  games  of  the  test.  Values  in  each  yield 
table  reflect  prior  decisions  on  the  frequency 
and  intensity  (table  1)  of  intermediate  cuttings. 

Computations  performed  by  the  subroutine 
follow  procedures  described  in  detail  elsewhere 
(Myers  1966,  Myers  and  Godsey  1968).  Average 
stand  diameter  (d.b.h.)  and  number  of  trees 
per  acre  just  before  initial  thinning  are  esti- 
mated from  measurements  of  stands  of  suitable 
ages  and  densities.  Number  of  trees  and  d.b.h. 
are  used  to  compute  initial  and  subsequent 
basal  areas  per  acre.  Total  cubic  volumes  per 
acre  are  computed  from  basal  area,  d.b.h.,  and 
tree  height  (table  2)  by  means  of  a  stand 
volume  equation  (statements  for  TOTO  and 
TOTT).  Total  cubic  volumes  are  multiplied 
by  factors  to  convert  them  to  merchantable 
cubic  volumes  (table  3)  and  board  feet  (table 
4). 

Increase  in  d.b.h.  due  to  initial  thinning  is 
computed  by  the  iterative  process  of  FORTRAN 
statements  6  to  12,  inclusive.  The  equation 
for  PDBHE  provides  an  estimate  of  diameter 
after  thinning  when  diameter  before  thinning 
and  the  percentage  of  trees  retained  are  known 
(Myers  1966).  Successive  percentages  are  tested 
until  d.b.h.  after  thinning,  number  of  trees, 
and  basal  area  agree  with  the  desired  thinning 
intensity  entered  as  THIN  by  INPUT1  and  as 
defined  in  table  1.  Equations  for  DBHP,  such 
as  statement  7,  express  that  part  of  table  1 
where  basal  area  is  less  than  66.3  square  feet. 

Table  1  gives  residual  basal  area  after  thin- 
ning for  various  average  stand  diameters.  The 
values  represent  one  possible  series  of  densities 
that  could  be  used  to  guide  successive  thinnings 
in  a  stand.  The  series  is  labeled  growing  stock 
level  80  to  indicate  that  reserve  basal  area  is 
80  square  feet  per  acre  when  d.b.h.  after  cutting 
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Table  1. — Basal  areas  after  intermediate  cutting  in  relation  to  average  stand  diameter 
Growing  stock  level  80  for  Black  Hills  ponderosa  pine 


Average 

Basal 

Average 

Basal 

Average 

Basal 

Average 

Basal 

stand  d.b.h. 

area 

stand  d.b.h. 

area 

stand  d.b.h. 

area 

stand  d.b.h. 

area 

after  cutting 

per 

after  cutting 

per 

after  cutting 

per 

after  cutting 

per 

( Inches ) 

acre 

( Inches ) 

acre 

(Inches) 

acre 

(Inches) 

acre 

Sq .   f t . 

Sq .    f t . 

Sq.  ft. 

Sq.  ft. 

2.0 

12.3 

4.0 

35.2 

6.0 

56.6 

8.0 

72.5 

2.1 

13.3 

4.1 

36.4 

6.1 

57.6 

8.1 

73.2 

2.2 

14.4 

4.2 

37.6 

6.2 

58.5 

8.2 

73.8 

2.3 

15.5 

4.3 

38.7 

6.3 

59.4 

8.3 

74.4 

2.4 

16.6 

4.4 

39.8 

6.4 

60.3 

8.4 

74.9 

2.5 

17.7 

4.5 

41.0 

6.5 

61.1 

8.5 

75.4 

2.6 

18.9 

4.6 

42.2 

6.6 

62.0 

8.6 

75.8 

2.7 

20.1 

4.7 

43.4 

6.7 

62.9 

8.7 

76.3 

2.8 

21.3 

4.8 

44.6 

6.8 

63.8 

8.8 

76.7 

2.9 

22.5 

4.9 

45.7 

6.9 

64.7 

8.9 

77.1 

3.0 

23.6 

5.0 

46.7 

7.0 

65.5 

9.0 

77.5 

3.1 

24.8 

5.1 

47.7 

7.1 

66.2 

9.1 

77.9 

3.2 

26.0 

5.2 

48.8 

7.2 

67.0 

9.2 

78.2 

3.3 

27.2 

5.3 

49.9 

7.3 

67.7 

9.3 

78.5 

3.4 

28.4 

5.4 

50.9 

7.4 

68.4 

9.4 

78.8 

3.5 

29.6 

5.5 

51.8 

7.5 

69.1 

9.5 

79.0 

3.6 

30.7 

5.6 

52.8 

7.6 

69.8 

9.6 

79.2 

3.7 

31.8 

5.7 

53.7 

7.7 

70.5 

9.7 

79.5 

3.8 

32.9 

5.8 

54.7 

7.8 

71.2 

9.8 

79.7 

3.9 

34.1 

5.9 

55.6 

7.9 

71.9 

9.9 

79.9 

10.0+  80.0 


Table  2. — Average  height  of  dominant  and  codominant  trees  at 
various  ages,  Black  Hills  ponderosa  pine 


Main 

Site  index 

class 

stand  age 

(Years) 

40 

1  50 

|  60 

70 

-  -  -  -  -  Feet 

10 

4.5 

4.5 

4.5 

4.5 

20 

9 

10 

12 

16 

30 

11 

16 

20 

26 

40 

17 

22 

28 

35 

50 

21 

28 

35 

43 

60 

26 

33 

41 

50 

70 

30 

38 

47 

56 

80 

34 

43 

52 

61 

90 

37 

47 

57  ~ 

66 

100 

40 

50 

60 

70 

110 

43 

53 

63 

74 

120 

45 

56 

66 

77 

130 

46 

59 

69 

80 

140 

48 

61 

71 

83 

150 

50 

63 

73 

86 

160 

51 

64 

75 

88 

170 

52 

65 

77 

90 

180 

53 

66 

78 

91 
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Table  3. — Factors  for  conversion  of  stand  volumes  in  total  cubic  feet  to 
merchantable  cubic  feet  per  acre,     Black  Hills  ponderosa  pine 


Average 

! 

Average 

Average 

stand 

Ratio 

of  | 

stand 

Ratio 

of  ! 

stand 

Ratio  of 

diameter 

merchantable  1 

diameter 

merchantable  j 

diameter 

merchantab le 

(Inches) 

to  total 

volume  : 

(Inches) 

to  total 

volume  | 

(Inches) 

to  total  volume 

5.0 

0.332 

8.1 

0.849 

11.9 

0.940 

5.1 

.355 

8.2 

.856 

12.  1 

.941 

5.2 

.  j  /  / 

8.3 

.  862 

12.4 

.942 

5.3 

.400 

8.4 

.868 

12.7 

.  943 

5.4 

.422 

8.5 

.872 

12.9 

.944 

5.5 

.444 

8.6 

.876 

13. 1 

.945 

5.6 

.465 

8.7 

.880 

13.3 

.946 

5.7 

.487 

8.8 

.884 

13.5 

.947 

5.8 

.508 

8.9 

.888 

13.  7 

.948 

5.9 

.530 

9.0 

.892 

13.9 

.949 

6.0 

.552 

9.1 

.896 

14.2 

.950 

6 .  1 

.575 

9.2 

.899 

14.4 

.951 

6.2 

.597 

9.3 

.  902 

14.  7 

.952 

6.3 

.618 

9.4 

.906 

14.9 

.953 

6.4 

.639 

9.5 

.910 

15.2 

.954 

6.5 

.659 

9.6 

.913 

15.4 

.955 

6.6 

.678 

9.7 

.916 

15.8 

.956 

6 . 7 

.  694 

9.8 

.920 

16.3 

.957 

6.8 

.710 

9.9 

.923 

16.8 

.  958 

6.9 

.725 

10.0 

.926 

17.3 

.959 

7.0 

.740 

10.1 

.928 

17.8 

.960 

7. 1 

.753 

10.2 

.930 

18.3 

.961 

7.2 

.766 

10.  3 

.931 

18.8 

.962 

7.3 

.778 

10.4 

.932 

19.3 

.963 

7.4 

.789 

10.5 

.933 

19.8 

.964 

7.5 

.799 

10.7 

.934 

20.3 

.965 

7.6 

.809 

10.9 

.935 

20.9 

.966 

7.7 

.818 

11. 1 

.936 

21.7 

.967 

7.8 

.826 

11.3 

.937 

22.5 

.  yoo 

7.9 

.834 

11.5 

.938 

23.3 

.969 

8.0 

.842 

11.7 

.939 

23.9 

.969 

1  To  4.0- inch  top 

in  trees 

6.0  inches  d.b.h. 

and  larger. 

Factor  for  an  unlisted  diameter  equals  factor  for  next  smaller  listed 
diameter.     For  example,   factor  for  15.6  inches  is  .955.^ 
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Table  4. — Factors  for  conversion  of  stand  volumes  in  total  cubic  feet  to  board  feet 
Scribner  rule  per  acre,1  Black  Hills  ponderosa  pine 


Average 

Ratio  of 

Average 

Ratio  of 

Average 

Ratio  of 

Average 

Ratio  of 

stand 

board  feet 

stand 

board  feet 

stand 

board  feet 

stand 

board  feet 

diameter 

to 

diameter 

to 

diameter 

to 

diameter 

to 

( Inches ) 

CllblC    f  £6t 

(Inches) 

cubic  feet 

( Inches ) 

cubic   f  S£t 

( Inche  s ) 

cubic  feet 

8.0 

0.  78 

11.9 

3. 

49 

15.  8 

4.71 

19.  7 

5 . 

42 

8.1 

.85 

12 

0 

3.56 

15.  9 

4.73 

19.8 

5 . 

44 

8.2 

.  92 

12 

1 

3. 

61 

16.0 

4.76 

19.  9 

5. 

45 

8.3 

.99 

12 

2 

3. 

65 

16.1 

4.78 

20.0 

5. 

46 

8.4 

1.06 

12 

3 

3.69 

16.2 

4.81 

20.1 

5. 

47 

8 . 5 

1. 13 

12 

4 

3. 

73 

16.3 

4.83 

20.  2 

5. 

48 

8.  6 

1.20 

12.5 

3. 

77 

16.4 

4.86 

20.  3 

5. 

50 

8.7 

1.  27 

12 

6 

3. 

80 

16.  5 

4.88 

20.4 

5. 

51 

8.8 

1.34 

12 

7 

3. 

84 

16.6 

4.90 

20.5 

5.52 

8.9 

1.41 

12 

8 

3. 

88 

16.7 

4.92 

20.6 

5. 

53 

9 . 0 

1 . 48 

12 

9 

3. 

91 

16.8 

4.  94 

20.7 

5 

54 

9 . 1 

1 . 55 

13 

0 

3. 

95 

16. 9 

4.96 

20.8 

5 , 

56 

9.  2 

1 . 62 

13 

1 

3. 

98 

17 . 0 

4.  98 

20.  9 

5 

57 

9.3 

1.68 

13 

2 

4. 

02 

17.1 

5.00 

21.0 

5. 

58 

9.4 

1.75 

13 

3 

4. 

05 

17.2 

5.02 

21.1 

5. 

59 

9.5 

1.82 

13 

4 

4. 

08 

17.3 

5.04 

21.2 

5. 

60 

9.6 

1.  89 

13 

5 

4. 

11 

17.4 

5.06 

21.3 

5. 

61 

9.  7 

1.  96 

13 

6 

4. 

14 

17.5 

5.08 

21.4 

5 

62 

9.8 

2.03 

13 

7 

4. 

17 

17.6 

5.10 

21.5 

5 

63 

9.9 

2.10 

13 

8 

4 

20 

17.7 

5.12 

21.6 

5 

64 

10.  0 

2.17 

13 

9 

4. 

23 

17.8 

5.13 

21.7 

5 

65 

10 . 1 

2 . 24 

14 

0 

4. 

25 

17.9 

5.15 

21.8 

5 

66 

10 . 2 

2 . 31 

14 

1 

4. 

28 

18.0 

5.17 

21.9 

5 

67 

10.3 

2.38 

14 

2 

4 

31 

18.1 

5.19 

22.0 

5 

68 

10.4 

2.45 

14 

3 

4. 

34 

18.2 

5.21 

22.1 

5 

69 

10 . 5 

9  ^9 
Z  .  JZ 

14 

4 

4. 

37 

18.  3 

5.22 

22.2 

5 

70 

10.  6 

Z  .  j  y 

14 

5 

4 

39 

18.4 

5.24 

22.3 

5 

71 

10.  7 

9  A1^ 
Z  .  DJ 

14 

6 

4 

42 

18.5 

5.26 

22.4 

5 

72 

10.8 

2.72 

14 

7 

4 

44 

18.6 

5.27 

22.5 

5 

73 

10.9 

2.79 

14 

8 

4 

47 

18.7 

5.29 

22.  6 

5 

74 

XX  .  U 

Z  .  oo 

14 

9 

4 

49 

18.8 

5.30 

22.7 

5 

75 

ii  l 

XX  .  X 

9 

z .  y  j 

15 

0 

4 

52 

18.  9 

5.32 

22.8 

5 

76 

11.2 

3.00 

15 

1 

4 

54 

19.0 

5.33 

22.9 

5 

77 

11.3 

3.07 

15 

2 

4 

56 

19 . 1 

j .  J  j 

9  ^  n 

Z  J  .  u 

5 

78 

11.4 

3.14 

15 

3 

4 

58 

19.2 

5.36 

23.1 

5 

79 

11.5 

3.21 

15 

4 

4 

61 

19.3 

5.37 

23.2 

5 

80 

11.6 

3.28 

15 

5 

4 

64 

19.4 

5.39 

23.3 

5 

81 

11.7 

3.35 

15 

6 

4 

66 

19.5 

5.40 

23.4 

5 

82 

11.8 

3.42 

15 

7 

4 

68 

19.6 

5.41 

23.5 

5 

83 

To  8-inch  top  in  trees  10.0  inches  d.b.h.  and  larger. 
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is  10  inches  or  larger.  Other  stocking  levels 
are  named  the  same  way.  For  example, 
stocking  level  100  means  that  reserve  basal 
area  will  be  100  square  feet  when  d.b.h.  after 
cutting  is  10  inches  or  larger.  Basal  areas 
for  level  100  and  for  diameters  smaller  than 
10  inches  are  obtained  by  multiplying  each 
basal  area  of  level  80  by  the  amount  100/80. 
Values  for  other  stocking  levels,  perhaps  from 
50  to  160,  are  computed  similarly. 

Periodic  increases  in  d.b.h.  due  to  tree 
growth  are  estimated  by  the  equation  for  DBHO 
in  the  loop  headed  DO  23.  The  equation  used 
in  the  example  is  for  Black  Hills  ponderosa 
pine,  and  for  a  projection  period  of  10  years. 
Equations  for  other  species  orprojectionperiods 
may  be  inserted  as  desired.  Intervals  between 
intermediate  cuttings  are  one  or  more  projection 
periods  long. 

It  is  often  desirable  to  make  simulations 
more  realistic  through  introduction  of  variability 
in  values  estimated  by  equations  or  contained 
in  tabulations.  For  example,  repeated  compu- 
tations of  DBHO  without  change  in  values  of 
the  independent  variables  will  always  give  the 
same  numerical  result.  In  reality,  actual  and 
estimated  values  frequently  differ.  A  way  of 
providing  variability  in  estimates  of  DBHO  is 
contained  in  the  program  segment  between 
statements  100  and  110.  Similar  statements 
could  be  written  for  other  variables. 

Variability  is  obtained  in  three  steps;  (1) 
generation  of  a  pseudorandom  number,  (2)  use 
of  this  number  as  an  independent  variable  to 
compute  the  value  of  a  residual  (range:  -0.3 
to  +0.3  inch),  and  (3)  addition  of  the  residual 
to  the  computed  value  of  DBHO.  The  pseudo- 
random number  generator,  statement  100,  is  of 
the  form: 

X.  =  AXi_1  +  C  (modulo  M) 

(Greenberger  1961).  Values  of  all  elements  of 
the  generator  are  specified  except  for  Xj_j  , 
which  is  read  in  as  variable  GNTR.  The  state- 
ment to  compute  RES  is  an  empirical  distribu- 
tion function  obtained  by  fitting  a  polynomial 
to  the  normally  distributed  residuals  of  the 
DBHO  equation  (Evans  et  al.  1967).  An  approx- 
imation to  the  normal  distribution  function 
may  also  be  used  (Burr  1967). 

Rethinnings  increase  d.b.h.  an  average  of 
0.4  inch  when  stand  densities  approximate 
growing  stock  levels  most  likely  to  be  goals 
of  timber  management. 


The  program  can  be  adapted  readily  for 
simulations  with  species  other  than  the  one 
used  for  the  test  problem.  Replacement  of 
tables  1  to  4,  several  statements  of  YIELD,  and 
species  designations  in  two  table  headings  are 
all  the  substitutions  necessary  if  the  same 
projection  procedure  is  used.  Statements  to 
be  replaced  are  the  equations  for  TOTO,  TOTT, 
DBHP,  PDBHE,  and  DBHO.  Projections  based 
on  other  procedures,  such  as  direct  estimation 
of  volume,  can  be  written  as  a  new  subroutine 
YIELD.  It  is  necessary  only  to  transfer  values 
of  CFMO(I),  BDFO(I),  CFMC(I),  and  BDFC(I) 
to  the  next  subroutine,  and  to  make  appropriate 
changes  in  the  READ  statements  in  INPUT1. 

Subroutine  ANVOL 

Subroutine  ANVOL  is  called  once  each  test 
to  compute  volume  per  acre  for  each  year  from 
initial  thinning  to  maximum  stand  age.  Volumes 
in  cubic  feet  and  board  feet  are  computed  by 
linear  interpolation  and  printed  in  a  composite 
table.  Stand  ages  cannot  exceed  179  years 
unless  dimensions  of  the  180-location  arrays  of 
acres  and  annual  volumes  are  increased. 

The  last  few  statements  of  the  routine 
expand  the  array  of  volumes  cut  to  assign  the 
volume  of  each  intermediate  cut  to  each  of  the 
years  before  the  next  cut  is  made.  These 
amounts  are  added  to  potential  postthinning 
volumes  in  subroutine  YEARS  to  compute 
volumes  per  acre  of  any  stands  not  given 
intermediate  cutting  because  they  are  older  than 
minimum  age  for  harvest. 

Subroutine  INPUT2 

INPUT2  is  called  once  each  game  to  enter 
numerical  values  of  variables  that  may  differ 
for  each  game  of  a  test.  Descriptive  data 
include  area  of  the  working  circle,  distribution 
of  area  by  age  classes,  nonstocked  area,  and 
number  of  acres  to  be  planted  annually. 
Various  costs  and  the  rate  at  which  they 
change  from  year  to  year  are  entered  into 
computer  memory.  One  to  ten  combinations 
of  price  limit,  allowable  cut,  and  minimum 
cutting  age  are  read  in  for  the  determination 
of  annual  cut  described  in  the  section  on  the 
data  deck. 

Values  read  by  INPUT2  that  do  not  appear 
in  other  tables  are  printed  as  a  record  of  game 
conditions  on  pages  headed  "alternatives  for 
this  game." 
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Subroutine  AREAS 

AREAS  is  called  once  each  game  to  compute 
volumes  and  area  distributions  at  the  end  of 
the  year  before  simulation  begins.  Acres  in 
each  1-year  age  class  are  expanded  to  obtain 
a  record  of  the  age  of  each  individual  acre. 
Total  area  (LAND)  cannot  exceed  1,000  acres 
unless  the  dimension  of  ACAGE(I)  is  increased. 
Age  of  the  oldest  acre  cannot  exceed  179  years 
unless  dimension  changes  listed  in  the  descrip- 
tion of  AN  VOL  are  made.  Acres  are  totaled 
by  10-year  age  classes.  A  table  of  initial  dis- 
tribution of  acres  by  1-year  and  10-year  age 
classes  is  printed. 

Growing  stock  volume  is  totaled  in  board 
feet  and  in  cubic  feet.  Volume  of  an  acre  will 
be  added  to  the  total  of  only  one  of  the  two 
volume  units.  The  unit  will  be  board  feet  if 
the  board-foot  volume  equals  or  exceeds  the 
value  of  variable  BFMRCH  read  by  INPUT1. 
No  volume  will  be  computed  for  the  acre  if 
stand  age  is  less  than  the  specified  minimum 
(AGMRCH). 

Values  of  volume  and  money  variables 
that  are  not  zero  at  start  of  simulation  are 
computed.  These  values  are  then  assigned 
space  in  one  of  two  2-dimensional  arrays  to 
preserve  them  for  printing  by  OUTPUT2. 


Subroutine  YEARS 

YEARS  is  called  each  year  of  each  game  to 
simulate  the  changes  in  volume  and  value 
produced  by  tree  growth,  cutting,  and  other 
events.  Changes  and  their  order  of  occurrence 
are  indicated  by  the  comment  statements  of 
the  source  program.  Several  items  are  described 
in  more  detail  below. 

A  specified  area  (IPLNT)  is  planted  each 
year  if  nonstocked  acres  exist.  Nonstocked 
acres  are  those  deforested  by  fire  or  other 
catastrophe,  and  do  not  include  harvested  acres 
that  will  restock  in  the  allotted  time.  Some 
or  all  harvested  acres  could  be  added  to  non- 
stocked  area  to  simulate  delays  or  failures  in 
natural  regeneration. 

Age  of  each  acre  destroyed  and  added  to 
nonstocked  area  is  selected  at  random  with  a 
pseudorandom  number  generator  of  the  form: 


Xi  =  AXi_1  +  C  (modulo  M) 


(Greenberger  1961).  All  values  are  preset  except 
for  X  j  _  j  which  is  designated  as  variable 
ANUL.  The  generator  has  a  periodicity  of 
128.  Values  of  ANUL  from  0  to  127  may  be 
selected  at  random  to  vary  the  pattern  of  loss. 

Annual  harvest  in  acres  equals  the  constant 
or  variable  allowable  cut  less  any  losses  of 
whole  acres.  Volume  and  value  of  shelterwood 
or  seed  trees,  if  any,  are  credited  to  the  year  of 
final  cut  FINL  years  after  the  regeneration 
cycle  starts.  The  volume  may  increase  or 
decrease  during  the  regeneration  period,  and 
may  be  left  unharvested.  Desired  results  are 
obtained  by  entry  of  appropriate  values  for 
GROW  and  SHELT. 

Subroutine  YEARS  was  written  to  contain 
a  series  of  dynamic  events  useful  for  many 
species  and  forest  regions.  The  computations 
can  easily  be  added  to  or  modified  to  meet 
local  needs  or  to  test  special  alternatives.  Un- 
wanted alternatives  in  the  program  need  not  be 
removed.  They  can  be  bypassed  by  entry  of 
appropriate  values  for  variables  not  needed. 


Subroutine  OUTPUT1 

Subroutine  OUTPUT1  is  called  after  the 
first  year  of  each  game  and  at  the  end  of  each 
decade.  Numbers  of  acres  by  1-year  and  10- 
year  stand  age  classes  are  printed.  The  tables 
correspond  to  that  printed  by  AREAS  just 
prior  to  start  of  a  game. 


Subroutine  ANUAL 

ANUAL  is  called  every  year  of  each  game 
to  compute  40  volume,  area,  or  money  totals 
and  to  store  them  for  later  use.  Each  total 
is  stored  in  one  of  two  2-dimensional  arrays. 
The  first  dimension  identifies  the  variable,  the 
second  the  year  of  a  game  to  which  the  value 
applies.  Numerical  value  of  each  year  subscript 
is  year  plus  one  so  year  zero  of  a  game  can 
be  included  in  the  array.  Array  values  are 
used  in  all  subsequent  subroutines  of  MAN AGD. 


Subroutine  OUTPUT2 

OUTPUT2  is  called  at  the  end  of  each  game 
to  print  the  results  of  each  year  of  the  game. 
Array  values  computed  and  stored  by  ANUAL 
are  printed  in  40  numbered  columns  that  extend 
across  four  pages  of  standard  Z-fold  paper. 
Entries  under  column  headings  are  printed  at 
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the  rate  of  40  lines,  or  years,  per  page.  Column 
headings  on  the  pages  produced  by  the  test 
problem  (appendix  2)  and  the  variable  lists  in 
the  source  program  of  ANUAL  (appendix  1) 
identify  the  variables  reported. 


Subroutine  WORTH 

Subroutine  WORTH  is  called  at  the  end  of 
each  game  to  discount  all  costs  incurred  and 
all  income  received.  Value  of  the  growing 
stock  at  the  end  of  the  simulation  period  is 
discounted  to  beginning  of  the  period.  The 
program  discounts  each  future  value  at  each 
of  20  compound  interest  rates.  Rates  range 
from  1.0  to  10.5  percent  at  intervals  of  0.5 
percent.  The  limits  and  interval  can  be  changed 
by  replacement  of  statements  for  CRATE(I) 
and  CRATE(K)  near  the  beginning  of  the 
subroutine. 

WORTH  prints  a  table  that  gives  the  present 
value  of  each  of  the  following  for  each  discount 
rate:  (1)  future  growing  stock,  (2)  all  incomes, 
(3)  sum  of  items  1  and  2,  (4)  all  costs,  and 
(5)  item  3  minus  the  sum  of  item  4  and  the 
value  of  the  growing  stock  at  beginning  of  the 
game.  Net  discounted  revenues  (present 
worths,  item  5)  may  be  plotted  over  discount 
rates  to  determine  the  internal  rate  of  return 
applicable  to  the  duration  and  conditions  of  the 
game. 


Subroutine  SUMRY 

Subroutine  SUMRY  may  be  called  at  the 
end  of  each  test  to  summarize  results  of  the 
games  of  the  test.  If  this  option  is  used, 
SUMRY  is  also  called  at  the  end  of  each  game 
to  store  specified  volume  and  money  values  in 
a  3-dimensional  array.  Values  stored  correspond 
to  the  columns  of  OUTPUT2  that  have  the 
column  numbers  entered  as  KOL(I)  by  INPUT1. 
Any  of  the  40  numbered  columns  of  OUTPUT2 
(appendix  2)  may  be  reproduced.  Not  more 
than  six  columns  may  be  summarized  during  one 
test  unless  the  dimensions  of  variables  KOL(I) 
and  SUMM(I,J,K)  are  increased.  Results  of  as 
many  as  10  games  may  be  summarized  at  one 
time. 

Results  of  the  games  of  a  test  are  printed 
together,  with  a  separate  page  for  each  variable 
selected  in  advance. 


Data  Deck 

Twenty-two  types  of  punch  cards,  listed 
below,  are  used  to  enter  initial  values  of  vari- 
ables into  computer  memory.  Most  cards  are 
not  optional  and  must  be  included  in  the  data 
deck  so  READ  statements  will  be  executed 
properly.  Four  types  are  optional  (6,  9,  10,  17) 
and  are  omitted  from  the  data  deck  if  the 
options  are  not  to  be  exercised. 

Data  cards  are  read  by  four  routines  in  the 
order  in  which  the  types  are  numbered.  Two 
types  are  read  once  by  the  main  program:  (1) 
card  type  1  or  BATCH  name,  and  (2)  card  type 
2  or  the  number  of  tests  to  be  performed  in  a 
batch.  These  identify  the  job  and  control  the 
number  of  times  the  rest  of  the  main  routine 
is  repeated. 

Card  types  3  to  10,  inclusive,  are  read  by 
INPUT1.  One  card  of  each  type  except  types 
6,  9,  and  10  must  be  read  once  each  test.  Card 
types  6  (1  card),  9  (15  cards),  and  10  (15  cards) 
are  omitted  from  the  data  deck  if  their  options 
are  not  to  be  used.  The  READ  statement  for 
card  type  6,  column  numbers  of  items  to  be 
summarized  by  SUMRY,  is  bypassed  when  zero 
is  punched  for  NKOLS  on  card  type  4.  Non- 
zero stumpage  prices  (BDPRI  and/ or  CFPRI) 
on  card  type  8  cause  the  corresponding  READ 
statements  for  variable  prices  of  card  types 
9  or  10  to  be  skipped. 

Subroutine  YIELD  reads  four  card  types 
once  each  test.  Types  11  to  14,  inclusive, 
contain  the  values  of  tables  1  to  4  of  this 
publication,  or  equivalent  information  applicable 
to  other  species  or  utilization  standards.  A 
type  consists  of  three  to  eight  punch  cards. 
Change  in  the  lowest  site  index  of  card  type 
12  (table  2)  or  in  the  minimum  diameter  of 
types  13  or  14  (tables  3  and  4)  will  require 
changes  in  the  statements  that  compute  array 
subscripts  and  probably  in  dimensions  of  the 
arrays. 

Card  types  15  to  22,  inclusive,  are  read  by 
INPUT2  once  each  game.  Each  type  consists 
of  one  card  except  for  optional  type  17,  which 
requires  10  punch  cards.  Statements  that  refer 
to  card  type  17,  variable  area  by  age  classes, 
are  bypassed  when  a  non-zero  value  is  punched 
for  KAREA  on  card  type  16. 

Card  types  18,  19,  and  20  contain  values  for 
the  price  control  procedure  of  Gould  and  O'Regan 
(1965).  The  number  of  acres  harvested  annually 
and  the  minimum  cutting  age  can  be  made  to 
vary  with  the  current  stumpage  price  of  saw 
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logs.  For  example,  in  the  second  game  of  the 
test  problem  (appendix  2),  5  acres  will  be  cut 
if  price  per  thousand  board  feet  does  not 
exceed  $12.00.  Seven  acres  will  be  cut  if  stump- 
age  price  is  $12.01  to  $15.00,  and  10  acres  will  be 
cut  if  price  exceeds  $15.00  but  is  less  than 
$99.00.  Sequence  of  harvest  is  from  oldest 
acre  to  youngest,  so  full  allowable  cut  will  be 


taken  only  if  sufficient  acres  above  minimum 
cutting  age  are  available.  If  price  control  is 
not  wanted,  entries  for  allowable  cut  in  columns 
1  to  4  of  card  type  19  and  for  cutting  age  in 
columns  1  to  8  of  card  type  20  are  the  desired 
constant  limits.  A  critical  price  greater  than 
the  largest  possible  price  (for  example,  $99.00) 
is  entered  in  columns  1  to  8  of  card  type  18. 


Order  and  Contents  of  the  Data  Deck 


Card  Read      Frequency    No.  of  Variable 

type    Optional        by  read  cards  name  Columns    Format  Description  of  variable 


1  NO  Main  Batch  1  BATCH (I)  1-24      3A8  Descriptive  name  to  identify  output  of 

one  pass  through  the  computer. 

2  NO  Main  Batch  1  NTSTS  1-4        14  Number  of  tests  in  the  batch,  each  with 

a  yield  table. 


3          NO  INPUTl        Test  1  DESCR(I)  1-40      5A8  Phrase  to  describe  conditions  of  one 

test;  to  identify  output. 


4         NO  INPUTl        Test  1          NGAME  1-4        14  Number  of  trials  (games)  to  be  operated 

in  one  test. 

NKOLS  5-8        14  Number  of  columns  of  OUT PUT 2  to  be 

printed  by  SUMRY. 

NOYRS  9-12      14  Number  of  years  simulated  in  each  game. 

Can  be  up  to  150,  but  will  usually  be 
less . 


INPUTl        Test  1         AGEO  1-8        F8.3        Stand  age  at  time  of  initial  thinning. 

First  age  given  in  the  yield  table. 

SITE  9-16      F8.3        Site  index.     Base  age  and  crown  class 

same  as  used  to  derive  growth  equations. 

DENO  17-24      F8.3        Number  of  trees  per  acre  at  age  AGEO. 

DBHO  25-32      F8.3        Average  diameter  breast  high  of  the 

stand  at  age  AGEO. 

ROTA  33-40      F8.3        Maximum  age  in  the  yield  table;   1  year 

more  than  the  maximum  age  expected  dur- 
ing simulations.     Cannot  exceed  180 
years . 

PRET  41-48      F8.3        Estimated  percentage  of  the  number  of 

trees  to  be  retained  in  initial  thinning 
at  age  AGEO.     Enter  as  a  percent,  e.g., 
35.0. 

DLEV  49-56      F8.3        Density  level  for  intermediate  cuts 

after  initial  thinning.     Based  on  table 
1  of  this  publication  and  procedure 
described  in  YIELD. 

CYCL  57-64      F8.3        Interval  between  intermediate  cuts. 

Equal  to  or  a  multiple  of  RINT. 

RXNT  65-72      F8.3        Number  of  years  for  which  a  growth  pro- 

jection is  made  by  the  equation  in  YIELD 

THIN  73-80      F8.3        Density  level  after  initial  thinning  at 

age  AGEO.     Based  on  table  1  and  proce- 
dure described  in  YIELD.     May  equal  DLEV. 
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Card  Read      Frequency    No.  of  Variable 

type    Optional        by  read  cards  name  Columns  Format 


Description  of  variable 


6  YES  INPUT 1  Test 


7         NO  INPUT1  Test 


NO  INPUTl  Test 


K0L(1) 

1-4 

14 

K0L(2) 

5-8 

14 

K0L(3) 

9-12 

14 

K0L(4) 

13-16 

14 

K0L(5) 

17-20 

14 

K0L(6) 

21-24 

14 

AGMRCH 

1-8 

F8, 

BFMRCH 

9-16 

F8, 

BFSALV 

17-24 

F8, 

COMCU 

25-32 

F8, 

COMBF 

33-40 

F8. 

BFPCT 

41-48 

F8, 

CFPCT 

49-56 

F8. 

GNTR 

57-64 

F8. 

BDPRI 

1-8 

F8 

CFPRI 

9-16 

F8 

Numbers  of  the  columns  of  0UTPUT2  to  be 
printed  by  SUMRY.     Maximum  number  of 
columns  is  6  and  must  agree  with  NKOLS 
of  card  type  4.     Column  numbers  are  from 
1  to  40,  as  given  in  the  column  headings 
of  the  printout  of  0UTPUT2  of  the  test 
problem. 

Minimum  stand  age  for  an  acre  to  be  in- 
cluded in  growing  stock  volume. 

Minimum  volume  in  M  bd.  ft.  for  an  acre 
to  be  included  in  board-foot  growing 
stock  volume. 

Minimum  volume  per  acre  in  M  bd.  ft.  for 
commercial  salvage  after  fire,  wind,  or 
other  loss. 

Minimum  cut  per  acre  in  merchantable 
cubic  feet  (table  3)  for  a  cut  to  be  of 
positive  commercial  value. 

Minimum  cut  per  acre  in  M  bd.  ft. 
(table  4)  for  a  cut  to  be  of  positive 
commercial  value. 

Ratio,  as  a  decimal,  of  board-foot 
stumpage  values  of  thinnings  to  board- 
foot  stumpage  values  of  harvests. 

Ratio,  as  a  decimal,  of  cubic-foot 
stumpage  values  of  thinnings  to  cubic- 
foot  stumpage  values  of  harvests. 

Any  number  between  0  and  1023  used  to 
generate  random  element  of  the  increase 
from  DBHT  to  DBHO.     Enter  number  larger 
than  1024  to  bypass  this  step. 

Stumpage  price  per  M  bd.  ft.  of  final 
harvest  if  price  is  constant  for  all 
years  of  a  game.     Enter  zero  if  variable 
prices  will  be  entered  with  card  type  10. 

Stumpage  price  per  100  cubic  feet  of 
final  harvest  if  price  is  constant  for 
all  years  of  a  game.     Enter  zero  if  vari- 
able prices  will  be  entered  with  card 
type  9. 


9  YES  INPUTl  Test 


PRICF(I)  1-80       10F8.3     Stumpage  price  per  100  cubic  feet  of 

harvest  for  each  of  150  years.  Used 
when  CFPRI  equals  zero. 


10         YES  INPUTl  Test 


15  PRIBD(I)  1-80       10F8.3     Stumpage  price  per  M  bd.   ft.  of  harvest 

for  each  of  150  years.  Used  when  BDPRI 
equals  zero. 


11  NO  YIELD  Test 


TABLl(K)  1-63      21F3.1    Basal  area  after  thinning  in  relation  to 

stand  diameter.  Values  of  table  1  copied 
on  punch  cards.     Used  with  DLEV  and  THIN. 


12  NO  YIELD  Test 


TABL2 (K,L) 


1-75 


25F3.1    Tree  heights  by  age  and  site  index  class. 

Values  of  table  2  copied  on  punch  cards. 
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Card  Read    Frequency    No.   of  Variable 

type    Optional          by  read  cards  name  Columns    Format  Description  of  variable 

13          NO  YIELD          Test  8  TABL3 (K)  1-72      24F3.3    Factors  for  conversion  of  total  cubic 

feet  to  merchantable  cubic  feet.  Values 
of  table  3  copied  on  punch  cards. 

1A         NO  YIELD          Test  6          TABLA(K)  1-78      26F3.2    Factors  for  conversion  of  total  cubic 

feet  to  board  feet.     Values  of  table  4 
copied  on  punch  cards. 

15  NO  INPUT 2        Game  1  GMNAM(I)  1-24      3A8  Descriptive  name  to  identify  each  game 

of  a  test. 

16  NO  INPUT 2        Game  1          LAND  1-4        14  Total  acres  in  simulated  working  circle. 

Maximum  is  1,000  acres. 

MOLD  5-8        14  Age  of  oldest  stand  in  the  working  cir- 

cle at  start  of  a  game.     Maximum  is  179 
years . 

NONSTK  9-12      I 4  Number  of  acres  non-stocked  at  start  of 

a  game.     Does  not  include  acres  harvest- 
ed the  year  before  simulation  begins  if 
regeneration  will  take  place  in  the 
allotted  time. 

KAREA  13-16      14  Number  of  acres  in  each  1-year  age  class 

when  there  is  equal  area  in  each  class 
except  for  NONSTK. 

IPLNT  17-20      14  Number  of  acres  of  NONSTK  regenerated 

annually  by  direct  seeding  or  planting 
at  a  cost  of  CPLT  per  acre. 

17  YES  INPUT2        Game  10  IACRE(I)  1-72       1814        Acres  in  each  1-year  age  class  from  0  to 

not  more  than  179.    Use  if  constant  area 
KAREA  is  not  wanted.     Include  NONSTK  in 
IACRE(l)  as  well  as  on  card  type  16. 

18  NO  INPUT2        Game  1  PRIDIV(I)  1-80      10F8.3    Limiting  prices  used  to  determine  annual 

cut  in  acres  and  minimum  cutting  age. 

19  NO  INPUT2        Game  1         MALCUT(I)  1-40       1014        Allowable  annual  cut  in  acres.  May 

vary  with  PRIDIV. 

20  NO  INPUT2        Game  1          FMRCHD(I)  1-80      10F8.3    Minimum  cutting  age.     May  vary  with 

PRIDIV. 

21  NO  INPUT 2        Game  1  SHELT  1-8        F8.3        Volume  of  shelterwood  in  M  bd.   ft.  Enter 

zero  if  shelterwood  or  seed  trees  are 
not  retained. 

RATE  9-16      F8.3        Rate  of  annual  increase  in  costs.  Enter 

zero  if  constant  costs  are  desired. 
Otherwise,  enter  percentage  as  a  decimal. 

CPLT  17-24      F8.3        Cost  of  regenerating  1  acre  by  seeding 

or  planting. 

CTHN  25-32      F8.3        Cost  per  acre  of  precommercial  thinning 

with  stand  conditions  as  specified  for 
the  simulation. 

CLOSS  33-40      F8.3        Cost  per  acre  of  cleanup  after  loss  due 

to  fire,  wind,  etc.,  when  volume  that 
can  be  salvaged  is  less  than  BFSALV. 


-11- 


Card  Read      Frequency    No.  of  Variable 

type    Optional        by  read  cards  name  Columns  Format 


Description  of  variable 


41-48  F8.3 


Total  per  acre  for  1  year  of  the  annual 
costs  that  can  be  assessed  by  area. 


F8.3 


Total  of  the  costs  that  can  be  assessed 
against  each  100  cubic  feet  harvested. 


57-64  F8.3 


Total  of  the  costs  that  can  be  assessed 
against  each  M  bd.  ft.  harvested. 


GROW  65-72      F8.3        Percentage  increase  or  decrease  in 

shelterwood  volume  during  the  regenera- 
tion period.     Enter  as  a  decimal. 
Enter  -1.0  if  shelterwood  or  seed  trees 
will  not  be  harvested. 


FINL  73-80      F8.3        Number  of  years  between  harvest  cut  and 

removal  of  shelterwood  or  seed  trees. 
Enter  zero  if  not  to  be  removed. 


22         NO  INPUT2      Game  1  DEFOR  1-8        F8.5        Percentage,  as  a  decimal,  of  the  area  of 

forest  lost  annually  to  fire,  wind,  etc. 

ANUL  9-16      F8.5        Any  number  between  0  and  127  used  to 

begin  generation  of  pseudorandom  numbers 
that  represent  ages  ot  stands  lost  to 
fire  or  other  agency. 


Description  of  Test  Problem 

The  test  problem  that  follows,  (detailed  in 
appendix  2)  demonstrates  most  computations 
possible  and  the  printed  results  obtained.  It 
may  be  used  to  verify  accuracy  of  punching  of 
source  decks  and  compatibility  of  the  program 
with  locally  available  compilers.  Growth  projec- 
tions and  volume  conversions  are  based  on 
relationships  applicable  to  Black  Hills  ponderosa 
pine  (Myers  1966).  Data  for  costs,  stump  age 
prices,  and  other  items  are  hypothetical.  Results 
of  the  simulations  are  therefore  examples  only, 
and  do  not  apply  to  any  real  forest  area. 

Assume  an  area  of  915  acres  of  managed 
stands  that  range  from  0  (just  harvested)  to 
129  years  old.  There  are  7  acres  in  each  1-year 
age  class,  plus  5  acres  of  old  burn  and  wind- 
throw  that  must  be  seeded  or  planted.  Annual 
losses  to  fire,  wind,  and  other  agencies  average 
0.04  percent  of  the  forested  area.  Site  index 
of  all  acres  is  60  feet. 

Stands  will  be  regenerated  by  two-cut 
shelterwood,  and  will  be  thinned  several  times 
during  a  rotation.  Shelterwood  volume  will 
average  4,000  board  feet  per  acre,  and  will 


increase  an  average  of  3  percent  of  initial 
volume  each  of  the  10  years  before  the  final 
cut.  Stands  will  be  thinned  for  the  first  time 
when  they  are  30  years  old.  At  this  age,  there 
will  be  1,000  trees  per  acre  that  average  4.5 
inches  in  diameter.  Initial  thinning  will  be  to 
level  120,  or  120/80  times  the  basal  areas  in 
table  1.  Rethinnings  at  20-year  intervals  will 
be  to  level  100,  or  100/80  times  tabulated  basal 
areas.  Stands  are  not  expected  to  ever  become 
150  years  old. 

Potential  prices  of  two  products  have  been 
estimated  for  each  of  the  next  30  years.  The 
stumpage  price  of  100  cubic  feet  of  roundwood 
from  mature  trees  or  from  thinnings  is  expected 
to  be  $2.50  throughout  the  period.  Price  of  a 
thousand  board  feet  of  saw  logs  is  expected  to 
vary  annually,  as  shown  in  column  28  of  the 
printout  of  annual  results  (appendix  2).  Saw 
logs  from  thinnings  will  sell  for  85  percent  of 
the  price  of  logs  from  harvest  cuts.  A 
minimum  commercial  cut  of  saw  logs  will  be 
3,000  board  feet  per  acre,  except  that  1,500  board 
feet  may  be  salvaged  from  an  acre  after  fire 
or  other  catastrophe.  Minimum  commercial 
cut  of  roundwood  will  be  400  cubic  feet. 
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Current  value  of  the  growing  stock  will  be 
computed  only  for  stands  at  least  40  years  old. 
Value  will  be  computed  for  cubic  volume  for 
acres  with  less  than  1,500  board  feet.  Other- 
wise, board-foot  volumes  will  be  used. 

Present  costs  of  various  operations  are  as 
follows: 

Costs  per  acre- 
Seeding— $30 

Precommercial  thinning— $25 
Cleanup  where  salvage  is  not  pos- 
sible—$25 
Annual  costs— $0.20 
Costs  assessed  against  volume  sold— 
Per  100  cubic  feet— $0.05 
Per  thousand  board  feet— $1.56 

These  costs  are  expected  to  increase  at  a  rate 
of  1  percent  annually.  Resources  are  available 
to  seed  1  acre  each  year. 

Two  possible  means  of  setting  the  allowable 
annual  cut  are  to  be  tested.  One  alternative 
is  to  harvest  7  acres  annually,  regardless  of 
price  fluctuations.  Stands  less  than  130  years 
old  will  not  be  cut.  A  second  possibility  is  to 
harvest:  (1)  5  acres  if  stumpage  price  per 
thousand  board  feet  is  $12  or  less,  (2)  7  acres 
if  the  price  is  $12.01  to  $15,  and  (3)  10  acres 
if  price  exceeds  $15  per  thousand.  Stands  less 
than  130  years  old  will  not  be  cut  except  that 
the  minimum  age  will  be  120  years  when 
stumpage  price  exceeds  $15.  Periodic  production 
in  board  feet  and  net  income  will  be  compared. 
Values  needed  to  obtain  present  worths  will 
be  computed. 

Data  cards  to  enter  these  values  into  compu- 
ter memory  must  contain  the  alphameric 
characters  given  in  the  following  list.  Decimal 
points  are  shown  for  numbers  in  F-format  to 
indicate  the  way  in  which  percentages  and 
money  values  are  entered.  Spaces  between 
numbers  do  not  correspond  to  the  blank  columns 
of  the  data  cards.  Card  types  9  and  17  are  not 
included  in  the  data  deck  because  the  options 
that  require  them  will  be  bypassed. 

Test  conditions  and  results  of  the  simula- 
tions are  printed  on  seven  types  of  pages 
(appendix  2).  The  first  two  types,  (1)  a  yield 
table  and  (2)  tables  of  volumes  per  acre  for 
each  year  of  stand  age,  appear  once  because 
one  test  was  run.  Four  types  of  pages  are 
printed  for  each  of  the  two  games.  The 
seventh  type  of  page  appears  once  at  the  end 
of  the  printout  to  summarize  specified  results 
of  the  two  games. 


The  two  sheets  of  "alternatives  for  this 
game"  show  the  values  used  in  the  simulations, 
including  the  different  allowable  cuts  and  cutting 
ages  tested. 

Distributions  of  acres  by  age  classes  (page 
type  four)  appear  on  two  sets  of  pages,  one  set 
for  each  game.  Pages  for  year  zero  show  7 
acres  in  each  of  129  1-year  age  classes.  Age 
class  zero  has  an  additional  5  acres  of  non- 
stocked  area.  Acreages  are  the  same  for  both 
games,  because  initial  distributions  were  the 
same.  Similar  pages  are  printed  at  the  end  of 
the  first  year  of  each  game  and  at  the  end  of 
each  decade.  For  brevity,  only  the  page  printed 
after  the  thirtieth  year  of  each  game  is  repro- 
duced in  appendix  2.  After  30  years  of  simula- 
tion, losses  and  direct  seeding  have  modified 
the  pattern  of  7-acre  units.  In  addition,  area 
distribution  of  the  second  game  has  been 
changed  by  the  variable  annual  cuts. 

The  fifth  type  of  page  is  a  set  of  four  pages 
for  each  game.  Values  in  many  of  the  40 
numbered  columns  differ  between  games. 
Board-foot  volumes  are  unequal  because  of 
variations  in  annual  cuts  of  mature  timber 
during  the  second  game.  This  caused  money 
values  to  differ  from  those  reported  for  the 
first  game. 

A  page  of  discounted  money  values,  the 
sixth  type  of  page,  is  printed  for  each  game. 
Rate  of  return  was  about  the  same  for  both 
games.  Both  operations  were  profitable.  In 
addition,  the  forest  would  probably  be  in  good 
condition  to  produce  other  products,  especially 
recreation. 

Last,  specified  values  from  each  game  were 
printed  together  for  convenience  in  interpreta- 
tion of  results.  Total  volume  in  board  feet  of 
all  cuts  plus  growing  stock  (column  10)  was 
higher  after  30  years  where  equal  areas  were 
cut  each  year.  Differences  were  never  great; 
variable  annual  cuts  produced  the  larger 
volume  after  20  years.  Total  net  worth  (column 
40)  was  greater  where  annual  cuts  varied  with 
price,  except  for  several  of  the  earlier  years. 

It  must  be  emphasized  that  results  of  these 
or  other  simulations  depend  on:  (1)  duration  of 
the  games,  (2)  values  entered  for  the  various 
variables,  (3)  assumptions  made,  and  (4)  degree 
to  which  the  model  represents  reality. 

The  above  information,  additional  data,  and 
knowledge  of  local  conditions  would  help  the 
forest  manager  decide  how  he  might  best  con- 
duct his  business.  Money  yields  would  en- 
courage the  manager  to  vary  annual  cuts  in 
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Data  Deck  for  Test  Problem 


Card 

type 


Contents  of  Cards 


1 
2 
3 
4 
5 
6 
7 
8 
10 
11 
12 
13 
14 
15 
16 
18 
19 
20 
21 
22 
15 
16 
18 
19 
20 
21 
22 


TEST  PROBLEM 
1 

MANAGED,  THINNED  AGE  30 
2  2  30 

30.  60.        1000.  4.5 

10  40 

40.  1.5  1.5  400. 

0.0  2.50 
(Column  28  of  printout  by  0UTPUT2) 
(Table  1  of  this  publication) 
(Table  2  of  this  publication) 
(Table  3  of  this  publication) 
(Table  4  of  this  publication) 
EQUAL  AREAS  CUT  ANNUALLY 


915 
99. 
7 

130. 
4.0 


129 
0 
0 

0 


01 


.0004  21. 
VARY  CUT  WITH  PRICE 


915 
12. 
5 

130. 
4.0 


129 
15. 

7 

130. 


.01 


5 
99. 
10 
120. 
30. 


7 
0 
0 
0 

25. 


150. 


3.0 


1 

0 
0 
0 

25. 


1 
0 
0 
0 

25. 


40. 


.85 


.20 


.20 


100. 


.05 


.05 


20. 


1.0  2222. 


0 
0 
0 

1.56 


0 
0 
0 

1.56 


10. 


.30 


.30 


120. 


0 
0 
0 

10.0 


0 
0 
0 

10.0 


0004  21. 
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response  to  changes  in  stumpage  price.  Highly 
variable  annual  cuts  and  equally  variable  net 
incomes  would  suggest  that  additional  simula- 
tions be  run  to  test  other  alternatives.  Cost 
of  computer  time  need  not  restrict  the  manager 
in  his  search  for  information.  The  test  problem 
was  compiled  and  run  in  26  seconds  of  central 
processor  time  and  9  seconds  of  input-output 
time. 
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APPENDIX  I 
Program  MANAGD 


DEFINITIONS   OF  VARIABLES 

ACAGEII)    =   ONE    ACRE,    WITH   AGE   OF   NUMERICAL   VALUE    AND    IDENTIFIED  BY 

SUBSCRIPT 
ACCST   =    ANNUAL   COST   PER  ACRE 
AGEO   =    INITIAL    AGE    IN   YIELD  TABLE 

AGMRCH   =   MINIMUM  AGE   FOR    STAND   TO   BE    INCLUDED    IN   GROW  I  NG  STOCK 
ANBDF(I)    =   M   BD.    FT.    PER   ACRE   AT   END   OF    EACH  YEAR 
ANCUV(I)    =   CU.FT.    STANDING   PER    ACRE    AT    END   OF    EACH  YEAR 
ANNE  T   =    ANNUAL    NET  INCOME 

ANUL   =  NUMBER   BETWEEN  0  AND   127  USED  TO  START-  GENERATION  OF 

PSEUDORANDOM  NUMBERS 
BASC   =    BASAL    AREA   REMOVED   PER  ACRE 
BASO   =   BASAL   AREA  PER  ACRE   BEFORE  THINNING 
BAST   =    BASAL    AREA   PER   ACRE    AFTER  THINNING 
BATCH ( I )    =   JOB  NAME 

BDFC(I)    =   M   BD.    FT.    REMOVED   PER  ACRE 
BDFO(I)    =   M  BD.    FT.   PER  ACRE   BEFORE  THINNING 
BDFT(I)    =   M   BD.    FT.    PER   ACRE    AFTER  THINNING 
BOPRI    =   CONSTANT   STUMPAGE    PRICE   PER   M   BD.  FT. 
BFCST   =   COSTS   PER   M   BD.    FT.  HARVESTED 

BFMRCH   =   MINIMUM   VOLUME    TO   BE    INCLUDED    IN   BD.  FT .    GROWING  STOCK 
BFPCT   =    PCT.    TO  CONVERT   BD.FT.    PRICE   FOR  THINNINGS 
BFSALV   =   MINIMUM   BD.FT.    FOR   COMMERCIAL  SALVAGE 
CFMC(I)    =   MERCHANTABLE   CU.FT.    REMOVED   PER  ACRE 

=   MERCHANTABLE   CU.FT.    PER   ACRE    BEFORE  THINNING 
=  MERCHANTABLE  CU.FT.    PER  ACRE   AFTER  THINNING 
PCT.    TO  CONVERT  CU.FT.    PRICE   FOR  THINNINGS 
CONSTANT   STUMPAGE    PRICE    PER    100  CU.FT. 
COST   OF   CLEANUP   OF   VOLUME   NOT  SALVAGED 
MINIMUM   COMMERCIAL   CUT    IN   BD.  FT. 
MINIMUM  COMMERCIAL   CUT    IN   CU.  FT. 
PLANTING     COST   PER  ACRE 


CFMO  (  I  ) 
CFMT  (  I  ) 
CFPCT  = 
CFPRI  = 
CLOSS  = 
COMBF  = 
COMCU  = 
CPLT  = 

CRATE!  I)   =    INTEREST  RATES   FOR  DISCOUNTING 

CSTAC  =   ANNUAL  COSTS   BASED  ON  AREA 

CSTVL   =    ANNUAL   COSTS   FOR   VOLUME  HARVESTED 

CTHN   =   COST   PER   ACRE   OF    PR ECOMME RC I AL  THINNING 

CUCST  =   COSTS   PER   100  CUBIC   FEET  HARVESTED 

CUTAGE   =   MINIMUM  CUTTING  AGE 

CYCL   =    INTERVAL    BETWEEN    INTERMEDIATE  CUTS 

DBHO   =    AVERAGE    STAND   D.B.H.    BEFORE  THINNING 

DBHT   =   AVERAGE    STAND   D.B.H.    AFTER  THINNING 

DEFOR   =    PERCENTAGE! AS   DECIMAL)    OF   NUMBER   OF   ACRES   LOST  ANNUALLY 

DENC   =    TREES   REMOVED   PER  ACRE 

DENO  =   TREES   PER  ACRE   BEFORE  THINNING 

DENT   =  TREES  PER  ACRE  AFTER  THINNING 

OESCR(I)    =  DESCRIPTION  OF  TEST  CONDITIONS 

DISCI  I)    =   DISCOUNTED  VALUE   OF   FUTURE  COSTS 

DISG(I)    =  DISCOUNTED  VALUE   OF  GROWING  STOCK 

DISI(I)    =   DISCOUNTED   VALUE   OF   FUTURE  INCOMES 

DLEV  =  GROWING   STOCK   LEVEL   FOR   SECOND  AND  SUBSEQUENT  THINNINGS 
FINL    =   YEARS    BETWEEN   HARVEST    AND  REMOVAL   OF  SHELTERWOOD 
FMRCHD ( I )   =  MINIMUM  CUTTING  AGE   BASED  ON  PRICE 
GMNAM(I)    =   NAME   OF   THE  GAME 

GNTR   =    PSEUOORANOOM   NUMBER   GENERATOR.    VALUE   0   TO  1023. 

GROW   =   GROWTH   RATE   OF  SHELTERWOOD 

GSVALB  =  DOLLAR  VALUE   OF   BD.FT.    GROWING  STOCK 

GSVALC  =   DOLLAR  VALUE  OF  CU.FT.   GROWING  STOCK 

GVLBF   =   GROWING   STOCK   VOLUME,    M   BD.  FT. 

GVLCU  =   GROWING   STOCK   VOLUME,  CU.FT. 

HTSO   =    TREE   HEIGHT   BEFORE  THINNING 

HTST   =   TREE   HEIGHT    AFTER  THINNING 

IACKE(I)    =   ACRES   OF    WORKING   CIRCLE    IN   EACH    1-YEAR   AGE   CLASS,  AT 

START   OF  GAME 
IALCUT   =   NUMBER   OF   ACRES   ALLOWABLE    ANNUAL  CUT 
IGAME   =   NUMBER   OF  GAME 

IPLNT   =   NUMBER   OF   NON-STOCKED   ACRES   REGENERATED  ANNUALLY 
ISUM(I)    =   TOTAL   ACRES   EACH   10-YR   AGE  CLASS 
ITEST   =   NUMBER   OF  TEST 

IVAR(I.J)    =   VARIABLES   TO   BE   PRINTED   BY  0UTPUT2 

IYEAR   =   YEAR   WITHIN   RUN   OF    A  GAME 

KAREA   =    EQUAL   AREA   OF   EACH    1-YEAR   AGE  CLASS 

KOL(I)    =   COLUMN   NUMBER ( FROM  0UTPUT2 )    PRINTED   BY  SUMRY 

KOUNT   =   COUNT   OF   ACRES   HARVESTED,    PLUS  ONE 

LAND   =    TOTAL   ACRES    IN   SIMULATED   WORKING  CIRCLE 

LAST   =   NUMBER   OF    LAST   ACRE  HARVESTED 

MALCUT(I)   =   ANNUAL   ALLOWABLE  CUT  BASED  ON  PRICE 

MOLD  =  AGE  OF  OLDEST  ACRE  IN  WORKING  CIRCLE  AT  START  OF  A  GAME 
IJACREII)   =  ACRES  OF  WORKING  CIRCLE   IN   EACH   1-YEAR   AGE  CLASS, 

DURING   A  GAME 
NGAMF   =  NUMBER  OF  GAMES   PER  TEST 

NKOLS   =    NUMBER   OF   COLUMNS   OF   0UTPUT2    TO   BE   PRINTED   BY  SUMRY 

NONSTK  =  NONSTOCKED  AREA 

NOYRS   =   NUMBER   OF   YEARS    IN   A  GAME 

NT  ST  S   =   NUMBER   OF   TESTS    IN  BATCH 

PRFT  =   PERCENTAGE  OF   TREES  RETAINED  AFTER   INITIAL  THINNING 
PREV(I)    =  PRESENT  VALUE   OF  GROWING  STOCK  AND  INCOMES 
PRIBD(I)    =  VARIABLE   STUMPAGE   PRICE  PER  M  BD.  FT. 
PRICF(I)    =    VARIABLE    STUMPAGE   PRICE   PER    100  CU.FT. 
PRIDIV(I)    =    PRICES   USED   TO    SET  POLICY 
PWTH(I)    =   PRESENT  WORTH 

RATE   =   RATE   OF   ANNUAL    INCREASE    IN  COSTS 
RETHV   =   ANNUAL   RETURN   FROM   FINAL  HARVEST 
RETRN  =    ANNUAL    INCOME   FROM  STUMPAGE 
RET  TH   =   ANNUAL   RETURN   FROM  THINNINGS 

R I  NT   =  NUMBER  OF  YEARS   FOR  WHICH  GROWTH  PROJECTION    IS  MADE 

ROTA   =   LONGEST   POSSIBLE   ROTATION    IN   YIELD  TABLE 

SCPLT    =    TOTAL   ANNUAL   PLANTING  COST 

SCTHN   =    SUM   OF    PRECOMMERC I AL    THINNING  COSTS 

SHELT   =   M   BD.    FT.    RETAINED   AS   SEED   TREES   OR  SHFLTERWOOD 

SHUD   =   VOLUME   OF    SHELTERWOOD   AT    FINAL  CUT 

SITE    =    SITE  INDEX 

SUMM(I,J,K)    =   ARRAY   FOR   PRINTING   BY  SUMRY 


TCOST   =    TOTAL   ANNUAL  COSTS 

THIN   =   GROWING   STOCK   LEVEL   FOR    INITIAL  THINNING 

TOTC    =   TOTAL   CUBIC   FEET   REMOVED   PER  ACRE 

rOTO   =   TOTAL   CUBIC   FEET   PER   ACRE    BEFORE  THINNING 

TOTT   =  TOTAL  CUBIC   FEET   PER  ACRE   AFTER  THINNING 

VAR(I.J)    =   VARIABLES    TO   BE   PRINTED   BY  0UTPUT2 

VBHV   =   BOARD-FOOT   VOLUME    FROM  HARVESTS 

VB TH   =   BOARD-FOOT   VOLUME   FROM  THINNING 

VCHV   =   CUBIC-FOOT   VOLUME   FROM  HARVESTS 

VCTH   =   CUBIC-FOOT   VOLUME   FROM  THINNING 

VLBF   =    VOLUME   HARVESTED,    M   BD.  FT. 

VLCU   =   VOLUME   HARVESTED,  CU.FT. 

YRLUS  =   NUMBER  OF   ACRES  LOST  ANNUALLY 

COMMON  ACCST, AGEO, AGMRCH , ANBDF ( 181) , ANCUVI 181) , ANNE T , BATCH ( 3 ), BDFC 
1 ( 180) ,BDFO( 180) , BFCST, CFMCI 180) ,CFMO( 180), CSTAC,CSTVL, CUCST, CUTAGE 
2,DBH0,DEN0,DESCR( 5) , DL EV , FMRCHD ( 10) , GMNAM I  3 ) , GSVALB , GSVALC , GVLBF , 
3GVLCU, I  ACRE (180) , I ALCUT , I G AME , I  SUM ( 18 ), ITEST, I VAR ( 26 ,  1 50 ) , I  YEAR , 
4K0L (6),LAND,LAST,MALCUT(10),M0LD,NACRE(180),NGAME,NKOLS,Nl , NONSTK, 
5N(JYRS,PRET, PRIBDI 150) ,PRICF(150),PRIDIV( 10) , RETRN, ROTA, RATE, GROW, 
6SITE,SUMM(6,25, 10), TCOST, SHE LT.YR LOS, FINL, CLOSS, CTHN, CPLT, CFPCT, 
7VAR( 14,150) , VLBF, VLCU, ACAGE ( 1000) , CYCL, KOUNT, DEFOR, ANUL, MIX, IPLNT, 
8RINT, THIN, BFMRCH, BFSALV, COMCU, COMBF, BFPCT, GNTR 

READ    (5,1)  (BATCH(I),I=1,3) 

1  FORMAT  (3A8) 
READ   (5,2)  NTSTS 

2  FORMAT  (14) 

DO  8    ITEST=1, NTSTS 
CALL  INPUT1 

PRINT  YIELD  TABLE   AND  COMPUTE  VOLUME   FOR  EACH  YEAR  OF  STAND  AGE 


OPERATE   SYSTEM    FOR    DESIRED  NUMBER   OF  GAMES 


DO  8  I  GAME= 1 »  NG  AME 
CALL    I NPUT2 


CREATE    ACRES    IN   EACH   AGE  CLASS 


OPERATE   SYSTEM   FOR   DESIRED   NUMBER   OF  YEARS 

DO   7    I  YE  AR= 1 , NOYRS 
CALL  YEARS 

PRINT    ACRES    IN   EACH   AGE   CLASS   FOR   FIRST   YEAR   AND   END  OF 
EACH  DECADE 

IF    (  I  YEAR    .LE.    1 )    GO   TO  3 
IF    (UK    .EQ.    10)    GO   TO  4 

GO   TO  6 

3  UK   =  1 
GO  TO  5 

4  UK   =  0 

5  CALL  0UTPUT1 

6  UK    =    UK   +  1 
CALL  ANUAL 

7  CONTINUE 

PRINT   VOLUMES   AND  VALUES   FOR   EACH  YEAR 

CALL  0UTPUT2 
CALL  WORTH 

SUMMARIZE   DESIRED  NUMBER  OF  COLUMNS  OF  0UTPUT2 

IF    (NKOLS    .LE.    0)    GO   TO  8 
CALL  SUMRY 

8  CONTINUE 
CALL  EXIT 

END 

SUBROUTINE  INPUT1 

COMMON   ACCST, AGEO, AGMRCH, ANBDF  (  181 ) , ANCUVI 181 )  , ANNE T , BATCH  (  3 )  ,  BDFC 
1 ( 180) ,BDFO( 180) , BFCST, CFMCI 180) ,CFMO( 180 ), CSTAC , CSTVL , CUCST,  CUTAGE 
2, DBHO, DENO, DESCR(5),DLEV,F MR CHD(10),GMNAM(3), GSVALB, GSVALC, GVLBF, 
3GVLCU,  I  ACRE ( 180) , IALCUT, IGAME, TSUMI 18) , ITEST, IV AR (26, 150), IYEAR, 
4K0L ( 6 ), LAND, LAS T, MALCUT (  10 ), MOLD, N ACRE!  180  ),NGA ME, NKOLS.Nl, NONSTK, 
5N0YR S, PRE T, PRIBDI 150) ,PRICF(150),PRIDIV(10) , RE T RN , ROT A , RATE , GROW , 
6SITE,SUMM(6,25, 10) , TCOST , SHFLT , YRLOS , F I NL , CLOS S , C THN , CPL T , CF PC T , 
7VARI  14, 150) , VLBF, VLCU, AC  AGE ( 1000)  ,CYCL, KOUNT, DEFOR, ANUL, MIX, IPLNT, 
8R I  NT, THIN, BFMRCH, BFSALV, COMCU, COMBF, 6FPCT, GNTR 

SET    INITIAL    VALUES   OF  ZERO 

DO    1  1=1,6 

1  KOLI I )    =  0 
DO  2  1=1,150 
PRIBDI  I  I    =  0.0 

2  PRICF  I  I  )   =  0.0 
DO   3  1=1,181 
ANBDF  I  I  )    =  0.0 

3  ANCUVI  I  )   =  0.0 
DO   4  1=1,6 

DO  4  J=l,25 
DO   4   K=l,  10 

4  SUMM( I , J,K)    =  0.0 


(EAD   VALUES   THAT   DO   NOT   CHANGE    WITHIN   A  TEST 
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=  1,5) 


READ   (5,5)    (DESCRI I), 

5  FORMAT  (5A8) 

READ    (5,6)  NGAME,NKOLS»NOYRS 

6  FORMAT  (2014) 

READ    (5,7)    AGEO, SITE, DENO.DBHO, ROTA, PRE T , DL E V, C YCL , R I  NT , THIN 

7  FORMAT  (10F8.3) 

IF    (NKOLS    .LE.    0)    GO  TO  8 
READ    (5,6)    ( KOL ( I ) , 1=1, NKOLS) 

8  READ   (5,7)    AGMRCH , BFMRCH , B FS AL V , C OMCU , COMBF , BF PCT , C FPC T , GNTR 
READ    (5,7)  BDPRI.CFPRI 

CREATE    A   SERIES   OF   CONSTANT   OR   VARIABLE  PRICES 

IF    (CFPRI    .NE.    0.0)    GO   TO  9 
READ    (5,7)    ( PRICFI I ) , 1=1, 150) 
GO   TO  11 

9  DO    10  1=1,150 

10  PRICFI I )   =  CFPRI 

11  IF  (BDPRI  -NE.  0.0)  GO  TO  12 
READ  (5,7)  ( PRIBDI I ) , 1=1, 150) 
GO   TO  14 

12  DO    13  1=1,150 

13  PRIBDI I )    =  BDPRI 

14  RETURN 
END 

SUBROUTINE  YIELD 

COMMON  ACCST, AGEO, AGMRCH, ANBOF I  181 ) ,ANCUV( 181) ,ANNET, BATCH!  3)  ,BDFC 
1(180),BDFO(180),BFCST,CFMC(180) ,CFMO( 1 80 ) , C STAC , C ST VL , C UCST , CUT  AGE 
2,DBH0,DEN0,0ESCR(5 ) , DLE V , FMRCHD I  1 0 ) , GMN AM  I  3 ) , GS V ALB , GS V ALC , GVL BF , 
3GVLCU, I  ACRE! 180) , I ALCUT , I  GAME , I  SUM  I  18) , ITEST, IVARI26, 150) , I  YEAR, 
4K0L ( 6 ), LAND, LAST, MALCUTI 10), MO LD,N ACRE! 180 ),NGAME, NKOLS, NL.NONSTK, 
5N0YRS, PRE T, PRIBDI 150) , PR  I C F ( 1 50 ) , PR  I D I V  I  10), RETRN, ROTA, RATE, GROW, 
6SITE,SUMM(6,25, 10 ) , TCOST , SHELT , YRLOS , F I NL , C LOS S , CTHN , C PLT , C FPC T , 
7VARI 14, 150) ,VLBF,VLCU,ACAGE( 1000)  , C YCL , KOUNT , DE FOR , ANUL , M I  X , I  PL  NT , 
8RINT, THIN, BFMRCH, BFSALV,COMCUiCOMBF,BFPCT, GNTR 

DIMENSION  TABL1 (81 ) ,TABL2(4,18) , T ABL3 I  192 ) , T ABL4 I  1 56 ) 

BDFT   =  0.0 

CFMT   =  0.0 

JBDFC   =  0 

J8DF0  =  0 

JBDFT   =  0 

JCFMC  =  0 

JCFMO  =  0 

JCFMT   =  0 

DO   1  1=1,180 

BDFCI I )  =0.0 

BDFOI I )    =  0.0 

CFMCI I )    =  0.0 

1  CFMOI I )    =  0.0 

READ  HEIGHTS, DENSITIES, VOLUME  CONVERSION  FACTORS 

READ    (5,2)  (TABL1IK),K=1,81) 

2  FORMAT  (21F3.1) 

READ    (5,3)    ( (TABL2(K,L),K=1,4),L=1,18) 

3  FORMAT  (25F3.1) 

READ    (5,4)  (TABL3(K),K=1,192) 

4  FORMAT  (24F3.3) 

READ    (5,5)    (TABL4IK) ,K=1, 156) 

5  FORMAT  (26F3.2) 
Nl   =  AGEO 

N  =  AGEO 

BASO   =   DENO   »   0.0054542    »   DBHO   *  DBHO 
OBTAIN  HTSO 
TEN  =  10.0 

ISITE   =    (SITE/TEN  -   3.0)    +  0.01 
I AGEO  =    (AGEO/TEN)    +  0.01 
HTSO  =   TABL2I ISITE, I AGEO) 

TOTO   =    (0.4047   »   BASO   •   HTSO)    +    (25.5970   »   DBHO)    -  191.6433 

COMPUTE   MERCH.    CU.FT.    IF   DBH    IS   AT   LEAST   5.0  INCHES 

IF    (DBHO    .LT.    5.0)    GO   TO  6 

I DBHO  =    ((DBHO  -   5.0)    »   10.0)    +  1.01 

XDBHO   =   TABL3I I DBHO ) 

CFMO(N)    =   TOTO   •  XDBHO 

COMPUTE   BD.FT.    IF   DBH   IS    AT   LEAST   8.0  INCHES 

IF    (DBHO   .LT.    8.0)    GO   TO  6 

JDBHO  =    ((DBHO  -   8.0)    *   10.0)    +  1.01 

YDBHO  =   T  ABL4 ( JDBHO ) 

BDFO(N)    =   TOTO   *  YDBHO 

COMPUTE   DBH   AFTER    INITIAL  THINNING 


DO   11  J=l,100 

PDBHE   =   0.95462»ALOG101DBHO)  - 

DBHE   =    10.0   •*  PDBHE 

IDBHE   =   DBHE   •    10.0   +  0.5 

DBHE   =  IDBHE 

DBHE   =  DBHE/TEN 

DENE   =   DENO   »  (PRET/100.0) 

BASE   =    (0.0054542   *   DBHE   *  DBH 

BREAK   =  49.9  *  THIN/80.0 

IF    (BASE    .GT.    BREAK)    GO  TO  7 

DI3HP   =    (80.0/THIN)    *  (0.08733 

GO   TO  8 

DBHP   =    (80.0/THIN)    *  (0.10938 

IUBHP   =   DBHP   •    10.0   +  0.5 

DBHP   =    I DBHP 

DBHP   =  D8HP/TEN 

IF    I DBHP  -  DBHE)  9,12,10 

PRET   =   PRET   +   PRET   •  0.02 

GO   TO  11 


0.10640*ALOG10IPRET) 


BASE)    -  0.17858 


10  PRET  =   PRET  -  PRET   •  0.02 

11  CONTINUE 

12  DBHT   =  DBHE 

COMPUTE   VALUES   AFTER   INITIAL  THINNING 

JDBHT   =    ( I DBHT   -   2.0)    »    10.0)    +  1.01 

SOFT   =   TABL1 I JDBHT ) 

BAST   =    (THIN/80.0)    «  SOFT 

ENTER   LOOP  FOR   ALL   REMAINING  COMPUTATIONS   AND  PRINTOUT 


DO   26  1=1,100 

DENT   =   BAST/IO. 0054542 

HTST   =  HTSO 

TOTT   =    (0.4047   »  BAST 

COMPUTE   MERCH.  CU.FT. 


♦   DBHT   *  DBHT) 
*   HTST)    +    (25.5970  *   DBHT)    -  191.6433 
IF   DBH   IS   AT   LEAST   5.0  INCHES 


IF    (DBHT   .LT.    5.0)    GO  TO  13 
I DBHT   =    I (DBHT   -   5.0)    •  10.0) 
XDBHT   =   TABL3 ( I DBHT ) 
CFMT   =   TOTT   *  XDBHT 


1.01 


COMPUTE   BD.FT.    IF   DBH  IS   AT   LEAST   8.0  INCHES 

IF    (DBHT    .LT.    8.0)    GO  TO  13 

JDBHT  =    ((DBHT   -   8.0)  *   10.0)    +  1.01 

YOBHT   =   T ABL4 ( J  DBHT ) 

BDFT   =   TOTT   *  YDBHT 

CHANGE   MODE   AND   ROUND  OFF   FOR  PRINTING 


13  JDENO 
JHTSO 
JTOTO 
JBASO 
JCFMO 
JBDFO 
JBDFO 
JDENT 
JHTST 
J  TOTT 
JDENC 
JCFMT 
CFMT 


JB  AST 
JBASC 
JTOTC 
JCFMC 


DENO 
HTSO 

TOTO  +  0.5 
BASO  +  0.5 
CFMOIN)    +  0.5 
( BDFO ( N ) /TEN )    +  0.5 
JBDFO  »  10 
DENT 
HTST 

TOTT  +  0.5 
JDENO  -  JDENT 
CFMT  +  0.5 
JCFMT 

IF    (JCFMT   .GT.    JCFMO)    JCFMO  =  JCFMT 
CFMOIN)    =  JCFMO 
JBDFT   =    I BDFT/TEN )    +  0.5 
JBDFT   =   JBDFT   *  10 
BDFT  =  JBDFT 
BDFT  =   BDFT   »  .001 

IF    (JBDFT   .GT.    JBDFO)    JBDFO  =  JBDFT 
BDFO(N)    =  JBDFO 
BDFO(N)    =   BDFO(N)    •  .001 
BAST   +  0.5 
JBASO  -  JBAST 
JTOTO  -   J  TOTT 
JCFMO  -  JCFMT 
IF    (JCFMC    .LE.    0)    JCFMC   =  0 
CFMC(N)    =  JCFMC 
JBDFC   =   JBDFO  -  JBDFT 
IF    (JBDFC    .LE.    0)    JBDFC   =  0 
BDFC(N)    =  JBDFC 
BDFC(N)    =   BDFC(N)    *  .001 
IF    ( I    . GE .    2)    GO   TO  19 

WRITE   HEADINGS   FOR   YIELD  TABLE 

WRITE    (6,14)  SITE.CYCL.DLEV 

14  FORMAT  (  1H1,////,28X,81HYIELDS  PER  ACRE  OF  MANAGED,  EVEN-AGED  STAN 
IDS  OF  PONDEROS A  PINE  IN  THE  BLACK  HILLS/1H  ,35X,10HSITE  INDEX, F3.0 
2,1H, ,F4.0,19H-YEAR  CUTTING  C YC LE , 1H , , 1 5H   DENSITY   LEVEL  ,F4.0,///) 

WRITE  (6,15) 

15  FORMAT  ( 1H0,25X,38HENTIRE  STAND  BEFORE  AND  AFTER  THINNING, 28X,26HP 
1ERI0DIC   CUT   AND  MORTALITY) 

WRITE  (6,16) 

16  FORMAT  ( 1H0,9X,5HSTAND, 1  OX , 5HB AS AL , 3 X , 7H AVER AGE , 2 X , 7H A VER AGE , 3 X , 5H 
1T0TAL,3X,9HMERCHANT-,3X,9HSAWTIMBER,9X,5HBASAL,4X,5HT0TAL,3X,9HMER 
2CHANT-,3X,9HSAWTIMBER) 

WRITE  (6,17) 

17  FORMAT  (1H  , 1  OX , 3H AGE , 4X , 5HTRE E S , 3X , 4H AR E A , 4X , 6HD . B . H . , 3X , 6HHE I GHT 
1,2X,6HV0LUME,2X,11HABLE  VOLUME , 4X , 6H VOLUME , 3X , 5HTRE E S , 3 X , 4H ARE  A , 3 X 
2, 6H VOLUME, 2X.11HABLE   VOLUME , 4X , 6HV0LUME ) 

WRITE  (6,18) 

18  FORMAT  (  1H  , 8X , 7H ( YE ARS ) , 3X , 3HN0 . , 3 X , 6HS0. F T . , 4X , 3H I N. , 6X , 3HFT .  ,  4X 
1,6HCU.FT.,5X,6HCU.FT.,8X,3HMBF,5X,3HN0.,3X,6HSQ.FT.,2X,6HCU.FT.,5X 
2,6HCU.FT.,8X,3HMBF) 

19  WRITE    (6,20)    AGEO , JDENO , JBASO , DBHO , JHTSO , JTOTO , C FMO ( N ), BDFO ( N ) 

20  FORMAT  ( 1H0 , 9X , F4. 0 , 4X , I  5 , 2X , I  4 , 5X , F 5 . 1 , 5X , I  3 , 4X , I  5 , 6X , F 5 . 0 , 6 X , F 6. 
13) 

IF    (AGEO. GE. ROTA)    GO  TO  27 

WRITE  (6,21 )  AGEO, JDENT, JBAST , DBHT, JHTST, J TOTT, CFMT, BDFT, JDENC, J  8 A 
ISC, JTOTC, CFMC IN) , BDFC I N ) 

21  FORMAT  I 1H  , 9X , F4 . 0 , 4X , I  5 , 2X , I  4, 5X , F5 . 1 , 5X , I  3 , 4X , I  5 , 6X , F5 . 0 , 6X , F6. 
13.4X,  I5..3X,  I3.5X,  I4,6X,F4.0,8X,F5.3) 

COMPUTE   VALUES   FOR   EACH   PERIOD.    THIN   AS  SPECIFIED 


IK  =  CYCL/RINT 
DO   2  3  L-ltlK 
AGEO  =   AGEO  +   R I  NT 
N  =  AGEO 

IF    (AGEO   .    GT   .    ROTA)    GO   TO  27 
DBHO  =    1.0097»DBHT   +  0.0096*SITE 
I DBHO  =   DBHO   »   10.0  +  0.5 
D8H0  =    I DBHO 
DBHO   =  DBHO/TEN 


1.5766*ALOG10(BAST) 


-   17  - 


ADD  RANDOM  ELEMENT  TO  PREDICTED  DBHO.IF  DESIRED 

IF    (GNTR    .GT.    1024.0)  GO  TO  110 

IDIV   =    (17.0   »   GNTR   +  3.01/1024.0 

NGNTR   =  GNTR 

GNTR  =    (17   *   NGNTR   +   3)    -    1024  »  IDIV 

IF   (GNTR   .GT.    1000.0)  GO  TO  100 

IF    (GNTR    .LT.    0.0)    GO  TO  100 
Al   =  GNTR/100.0 


A2 


Al 


0.0523 


A2 


IRES 
IRES 


0.5 
0.5 


ADJ/TEN 


RFS  =  0.9565 
1*A2  -  3.3009 
IRES  =  RES 
IF    (RES  .LT. 
IF    (RES   . GT. 
ADJ   =  IRES 
DBHO   =  DBHO 
110   DENO  =  DENT 

8AS0  =   DENO   *    (0.0054542  » 
ISITE  =    (SITE/TEN  -  3.0)  + 
IAGEO  =   (AGEO/TEN)    +  0.01 
HTSO   =   TABL2I ISITE, IAGEO) 
TOTO  =   0.4047   »   BASO   »  HTSO 


COMPUTE   MERCH.    CU.FT.    IF   DBH    IS   AT  LEAST   5.0  INCHES 

IF    (DBHO   .LT.   5.0)    GO  TO  22 

I DBHO  =    ((DBHO  -   5.0)    *    10.0)    +  1.01 

XDBHO   =   TABL3 ( IDRHO) 

CFMO(N)    =   TOTO   •  XDBHO 

COMPUTE   BD.FT.    IF   DBH    IS    AT   LEAST   8.0  INCHES 

IF    (DBHO    .LT.    8.0)    GO   TO  22 
JDBHO  =    ((DBHO  -   8.0)    *    10.0)    +  1.01 
YDBHO  =   TABL4 ( JDBHO ) 
BDFO(N)    =   TOTO   »  YOBHO 
22    IF    (L    .EC.    IK)    GO   TO  24 

WRITE   VALUES   FOR   END   OF   PERIOD    IF   THINNING  NOT  DUE 


A2   +  0.00084 


DENO 
HTSO 

BASO  +  0.5 
TOTO  +  0.5 
CFMO(N)    +  0.? 
=  JCFMO 
(  BDFO ( N ) /TEN  I 
JBDFO  »  10 
=  JBDFO 

=   BDFO(N)    *  .001 

AGEO,KDENO,KBASO, DBHO, KHTSO.KTOTO, CFMO ( N ) , BDFO ( M) 


0.5 


KQENO  = 
KHTSO  = 
KBASO  = 
KTOTO  = 
JCFMO  = 
CFKOIN) 
JBDFO  = 
JBDFO  = 
BDFO ( N ) 
BDFO(N) 
WRITE  (6,20) 
DBHT  =  DBHO 
BAST  =  BASO 
23  CONTINUE 


INCREASE   DBH   AS   RESULT   OF   THINNING   AND   COMPUTE  POST-THINNING 
VALUES 

24  DBHT   =    DBHO   +  0.4 

IF    (DBHT    . GE .    10.0)    GO   TO  25 

JOBHT   =    ((DBHT   -   2.0)    *    10.0)    +  1.01 

SQFT   =   TABL 1 ( JDBHT ) 

BAST   =    (DLEV/80.0)    *  SQFT 

GO  TO  26 

25  BAST   =  DLEV 

26  CONTINUE 

27  RETURN 
END 

SUBROUTINE  ANVOL 

COMMON   ACCST, AGEO, AGMRCH, ANBDFI 1 8 1 1 , ANCUV ( 1 8 1 1 , ANNET , B ATCH ( 3 1 , BDFC 
1 ( ieO) ,BDFO( 180) ,BFCST,CFMC( 180) ,CFMO( 1 8 0 ) , C STAC , C STVL , CUC ST , CUT AGE 
2,DBHO,OENO,DESCR(5) , OLE V , FMRCHD ( 10) ,GMNAM(3) , GS VALB , GSVALC , GVLBF , 
3GVLCU, I  ACRE (180) , IALCUT, I  GAME, ISUMf 18),ITEST,IVAR(26,150) , I  YEAR, 
4K0LI6) , LAND, LAST, MALCUTI 10) .MOLD, NACRE (180) , NGAME , NKOLS , Nl , NONSTK , 
5N0YRS,PRET,PRIBD( 150),PRICF(150),PRIDIV( 10) ,RETRN, ROTA, RATE, GROW, 
6SITE,SUMM(6,25, 10 ) , TCOST , SHELT , YRLOS , F I NL , C LOSS , C THN , C P LT , CF PC T , 
7VAR(14,150),VLBF,VLCU,ACAGE(1000),CYCL,K0UNT,DEF0R,ANUL,MIX,IPLNT, 
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF,BFPCT, GNTR 

IROT  =  ROTA 

INT   =   R I  NT 

.NVUL   =    (  (  IROT   -  Nl  )  /INT)    +  1 
K  =   NVOL   -  1 


INTERPOLATE    BETWEEN   VALUES   FROM   YIELD  TABLE 


'DO    1  L=1,K 
DO    1  J=1,INT 

NN  =   J   +   Nl   +    ( L   -   1 )    «  INT 

RJ  =  J  -  1 

N   =   Nl   +    (L  -   1)    »  INT 

ANCUV (NN)  =CFMO(N)-CFMC(N)+(RJ/RINT)( 
ANBDFI NN)    = BDFO ( N 1 -BDFC ( N 1  +  ( R J / R I  NT  1 < 
1  CONTINUE 


(CFMO(N+INT)-CFMO(N)+CFMC(N) ) 
(BDFO(N+INT)-BDFO(N)+BDFC(N) ) 


WRITE    TABLE  HEADINGS 

WRITE    (6,2)    SITE, CYCL, THIN, DLEV 

2  FORMAT  ( 1H1,////,41X,51HGR0WING  STOCK  OF  MANAGED  BLACK  HILLS  PONDE 
1R0SA  PINE/1H  ,47X,10HSITE  INDEX, F3.0,1H,,F4.0,19H-YEAR  CUTTING  CYC 
2LE/1H    ,53X, 14HDENSITY   L E VEL- , F4 . 0 , IX , 3H AND , F4 . 0 ) 

WRITE  (6,3) 

3  FORMAT    ( 1H0,43X,44HV0LUMES   PRESENT   PER   ACRE   AT   END   OF   EACH  YEAR) 
WRITE  (6,4) 

4  FORMAT  ( 1H  ,54X,23HMERCHANTABLE  CUBIC  FEET/1H0,64X,4HYEAR/1H  ,14X, 
16HDECADE,9X,1H0,9X, 1H1.9X, 1H2 , 9X , 1H3 , 9X , 1H4 , 9X , 1H 5 , 9X , 1H6 , 9X , 1H7 , 9 
2X, 1H8.9X, 1H9,//  ) 


WRITE   CUBIC    FEET   PER   ACRE   FOR   EACH  YEAR 

WRITE    (6,5)    K, ( ANCUVI NN 1 ,NN=1 , 10) 
i   FORMAT    (1H    ,  120, F13. 1 ,9F10. 1 ) 
IJ   =    IROT/10  -  1 
DO   6   J  =  l ,  I J 
NN=10*J+1 

WRITE  (6,5)  J, ANCUV (NN) , ANCUV (NN+1) , ANCUV (NN  +  2 )  , ANCU V ( NN+3 ),  ANCUV  ( 
1NN  +  4)  , ANCUV (NN  +  5)  , ANCUV ( NN  +  6 ) , ANCUV ( NN  +  7 ) , ANCUV (NN  +  8) , ANCUV  (  NN  +  9  ) 

,  CONTINUE 
J   =  IROT/10 

ANCUVI IROT+l )    =  CFMO(IROT) 
WRITE    (6,5)    J , ANCUVI IROT+l ) 

WRITE    BD.FT.    PER   ACRE   FOR    EACH  YEAR 

WRITE  (6,7) 

'   FORMAT    I 1H0,///,55X,23HTH0USANDS   OF   BOARD  FEET,//) 

WRITE    (6,8)    K, ( ANBDF(NN) ,NN=1, 10) 
I    FORMAT    ( 1H  ,I20,F13.3,9F10.3) 

DO   9   J  =  1,1  J 

NN  =    10   »   J   +  1 

WRITE  (6,8)  J,ANBDF(NN)  ,ANBDF(NN  +  1)  ,ANB0F(NN  +  2)  ,ANBDF(NN  +  3), ANBDFI 
l.NN  +  4)  , ANBDF ( NN  +  5 ) ,ANBDF(NN+6) ,ANBDF(NN  +  7) ,ANBDF(NN  +  8) , ANBDF  (  NN  +  9  ) 

I  CONTINUE 
J   =  IROT/10 

ANBDFI IROT+l)    =   BDFOI IROT) 
WRITE    (6,8)    J, ANBDFI IROT+l  ) 

PROVIDE   FOR   ANY   ACRES  BEYOND  CUTAGE   LEFT  UNTHI NNED  FOR  A  FEW  YEARS 

JCYCL   =  CYCL 

IMI    =   JCYCL   -  1 

DO   10    I=N1, IROT, JCYCL 

DO   10  J=1,IMI 

NX  =   I   +  J 


180)   GO  TO  11 
BDFC  I  I  ) 
CFMCI  I  1 


IF    (NX  .GE. 

BDFC ( NX ) 

10  CFMC(NX) 

11  RETURN 
END 

SUBROUTINE  INPUT2 

COMMON  ACCST, AGEO, AGMRCH, ANBDFI 181) , ANCUV  I  181 ) , ANNET , BATCH  I  3  ),  BDFC 
1(180), BDFOI180) ,BFCST, CFMCI 180) ,CFMO( 180  ), CST AC , C STVL , CUC ST ,  CUTAGE 
2,DBH0,DEN0,DESCR( 5) , DLE V , F MRCHD ( 1 0 1 , GMN AM ( 3 1 , G S V A  LB , GS V ALC , GVLBF , 
3GVLCU, I  ACRE! 180) , IALCUT, I G AME , I  SUM ( 1 8 ) , I T EST , I VAR ( 26 ,  150) ,  I  YEAR, 
4K0LI6) , LAND, LAST, MALCUTI  10), MOLD, NACRE  I  180)  , NG AME , NKOLS , N 1 ,  NONSTK , 
5N0YRS, PRET, PRIBDI 150),PRICF(150),PRIDIV( 10  1  , RE TRN  ,  ROT A , RATE ,  GROW  , 
6SITE,SUMM(6, 25, 10), TCOST, SHELT, YRLOS, FINL,CLOSS,CTHN,CPLT,CFPCT, 
7VAR(14,150)  ,VLBF, VLCU, AC  AGE ( 1000) , C YCL , KOUNT , DE FOR , ANUL , M I X ,  I  PL  NT , 
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF, BFPCT.GNTR 

SET    INITIAL   VALUES   OF  ZERO 

CSTAC   =  0.0 
CSTVL  =  0.0 
IALCUT   =  0 
KOUNT  =  1 
LAST  =  0 
RETRN  =  0.0 
YRLOS  =  0.0 
DO   1  1=1,180 
I  ACRE ( I  )    =  0 

1  NACRE ( I )    =  0 
DO  2  1=1,10 
FMRCHD ( I )    =  0.0 
MALCUT(I)    =  0 

2  PRIDIVI I )    =  0.0 
DO  3  1=1,26 

DO  3  J=l,150 

3  IVARI  I  ,  J)    =  0 
DO  4  1=1,14 
DO  4  J=l,150 

4  VAR( I , J  )   =  0.0 

READ   VALUES   THAT   DO   NOT   CHANGE   WITHIN   A  GAME 

READ    (5,5)    (GMNAMI I ) , 1=1,3) 

5  FORMAT  (3A8) 

READ    (5,6)    LAND, MOLD, NONSTK, KAREA, I PLNT 

6  FORMAT  (2014) 

IF    (KAREA   .EO.    0)    GO   TO  8 

READ    IN   EQUAL   AREAS    IN   ALL    AGE  CLASSES 

NDX   =   MOLD   +  1 
DO   7    1=1, NDX 

7  IACRE(I)   =  KAREA 

ADJUST   NUMBER   OF   ACRES    IN   OLDEST   CLASS    IF   TOTAL    AREA   NOT  MULTIPLE 

OF  KAREA 

KDIFF   =   LAND  -   NDX   »   KAREA  -  NONSTK 
I  ACRE ( NDX )    =    I  ACRE  I NOX )    +  KDIFF 
IACRE(l)   =   IACRE(l)    +  NONSTK 
GO   TO  10 

READ    IN   UNEQUAL    AREAS    IN    AGE  CLASSES 


8  READ    (5,9)    ( I  ACRE ( I  1 ,1  =  1,180) 

9  FORMAT  (1514) 


READ   LIMITATIONS   ON  CUT 


10  READ    (5,11)    (PRIDIVI I  I 

11  FORMAT  (10F8.3) 


18  - 


READ    (5,9)    (MALCUTI I ) , 1=1, 10) 
READ    (5,11)    ( FMRCHOI I ) , 1=1, 10) 

READ  (5,11)  SHE LT,R ATE, CPLT,CTHN,CLOSS, ACCST, CUCST, BFCST, GROW, F I NL 
READ    (5,12)    DEFOR, ANUL 

12  FORMAT  (2F8.5) 
DEFGR1    =   DEFOR   *  100.0 

PRINT   CONDITIONS   OF  SIMULATIONS 

WRITE  (6,13) 

13  FORMAT    (  1H1  ,//,46X,26HALTERNATIVES   FOR   THIS  GAME) 
WRITE    (6,14)    ( BATCH! I ), 1=1, 3) 

14  FORMAT    (1H    ,45X,7HRATCH  ,3A8) 
WRITF    (6,15)  ITEST 

15  FORMAT    (1H    , 45X , 4HTE S T , I  4 ) 
WRITE    (6,16)    ( GMNAM (I), 1=1, 3) 

16  FORMAT    ( 1H    ,45X,6HGAME  ,3A8) 
WRITE    (6,17)    (DESCRl I  )  , 1  =  1  ,5) 

17  FORMAT    (1H    , 4 5 X , 5  A 8 , / / / / ) 
WRITE    (6,18)  NOYRS 

18  FORMAT    (1H    , 45X , 24HNUMBER  OF   YEARS   PER   GAME, 14,//// ) 
WRITF    (6,19)    ( PRIDIV( I ) , 1=1, 10) 

19  FORMAT    (1H    .15HCRITICAL    PR  I C E S , 1 2 X , 1 0F9 . 2 ) 
WRITE    (6,201    (MALCUTI II , 1=1, 10) 

20  FORMAT    ( 1H   , 13HALL0WABLE   C UT , 1 1 X , 10  I  9 ) 
WRITF    (6,21)    (FMRCHDI I) , 1=1, 10) 

?1    FORMAT    ( 1H    , 19HM I N I  MUM  CUTTING  AGE , 5X , 1  OF  9. 0 , / / / / ) 
WRITE    (6,22)  LAND 

22  FORMAT    (1H    , 23HACRES    IN  WORKING  CIRCLE, 13X, I4,25X,27HC0STS   IN  FIRS 
IT   YEAR   OF  GAME) 

WRITE    (6,23)  ACCST 

23  FORMAT    (1H    ,69X,17HPER   ACRE    (ANNUAL), 8X,F9. 2) 
WRITE    (6,24)  CUCST 

24  FORMAT    (1H    , 38HM I N I  MUM  VALUES   FOR    INCLUSION    IN   TOTAL S , 3 1 X , 2 5HPER  1 
100   CU.    FT.    HARVESTED, F9. 2) 

WRITE    (6,25)  AGMRCH.BFCST 

25  FORMAT    ( 1H    ,4X,22HAGE,    FOR   GROWING   STOCK , 1 1 X , F3 . 0 , 2 9X , 1 3HPER   M  BD. 
1    FT. , 12X.F9.2) 

WRITE    (6,26)    BFMRCH, CTHN 

26  FORMAT    (1H    ,4X,28HM   BD.    FT.,    FOR   GROWING   STOCK , 5X , F5 . 1 , 27X , 1 3HTH I N 
1   ONE  ACRE,12X,F9.2) 

WRITE    (6,27)  COMCU.CPLT 

27  FORMAT    ( 1H    ,4X,27HCU.    FT.,    FOR   COMMERCIAL   CUT , 5X , F4. 0 , 29X , 1 4HPL A  NT 
1   ONE  ACRE,11X,F9.2) 

WRITE    (6,28)  COMBF.CLOSS 

28  FORMAT    (1H    ,4X,29HM   BD.    FT.,    FOR   COMMERCIAL   CUT , 4X , F 5 . 1 , 27X , 19HC LE 
1ANUP   OF   ONE  ACRE,6X,F9.2) 

WRITE    (6,29)    BF  S ALV , RATE 

29  FURMAT    ( 1H   .4X.22HM   BD.    FT.,    FOR    S AL V AGE , 1 1 X , F 5 . 1 , 23X , 2 5HR ATE   OF  I 
INCREASE    IN  COS T S , 4X , F9 . 2 , / / ) 

WRITE    (6,30)  IPLNT 

30  FORMAT    (1H   , 22HACRES   PLANTED   ANNUALL Y , 1 4X , 1 4 , 2 5X , 3 5HREL AT  I VE  VALUE 
1    OF    INTERMEDIATE  CUTS) 

WRITE    (6,31)  DEFORl.CFPCT 

31  FORMAT    (1H    , 30HPERCENT   OF   ACRES   LOST    ANNU ALL Y , 6X , F8 . 3 , 2 5X , 2 3HSTUMP 
1  AGE    PRICE,    CD.  FT.,2X,F9.2) 

WRITE   (6,32)  SHELT.BFPCT 

32  FORMAT    ( 1H    ,24HM  BD.    FT.    IN   SHELT ERWOOD , 1 3X , F5 . 1 , 27X , 2 3HSTUMP AGE  P 
lRICh,    BD.  FT.,2X,F9.2) 

WRITE    (6,33)  FINL 

33  FORMAT    (1H    , 19HREGENERAT ION   PER IOD, 18X , F 5 . 1 , // ) 
WRITE    (6,34)  ANUL 

34  FORMAT    (1H    , 29HPSEUD0RAND0M  NUMBER  GENER ATOR , 5 X , F 8 . 1 ) 
WRITE    (6,35)  GNTR 

35  FORMAT    (1H  ,34X,F8.1) 
RETURN 

END 

SUBROUTINE  AREAS 

COMMON   ACCST, AGEO, AGMRCH, AN BDF (181),ANCUV(181),ANNET,BATCH(3),BDFC 
1( 180) ,BDFO( 180) ,BFCST,CFMC( 180) ,CFMO( 1 80 ) , C S T AC , C ST VL , CUC ST , CUT AGE 
2, DBHO.OENO, DESCRl 5 ) , DLE V, FMRCHD ( 10) ,GMNAM( 3) , G S VA LB , GS V ALC , G VLBF , 
3GVLCU,IACRE(180),IALCUT,IGAME,ISUM(18), ITEST, IVAR(26, 150), IY EAR, 
4K0L (6) .LAND, LAST, MALCUT( 1 0 ) , MOL D, NACRE ( 1 80 ) , NG AME , NK OL S , N 1 , NONS TK , 
5N0YRS,PRET,PRIBD( 150) ,PRICF(150),PRIDIV(10)  , RE TRN , ROT  A , R ATE , GROW , 
6SITE,SUMM(6,25, 10) , TCOS T , SHE LT , YRLOS , F I NL , CLOS S , CTHN , CPLT , CFPC T , 
7VAR( 14, 150) ,VLBF,VLCU, AC  AGE ( 1 000 ) , C YCL , KOUNT , 0 EFOR , ANUL , MI  X , I PLNT , 
8RINT, THIN, BFMRCH, BFSALV, CO MCU,C0MBF,8FPCT, GNTR 

GSVALB  =  0.0 

GSVALC   =  0.0 

GVLBF   =  0.0 

GVLCU  =  0.0 

DO    1  1=1,18 

1  ISUMt I )    =  0 
DO   2  1=1,1000 

2  ACAGE ( I )    =  0.0 

CONVERT   ACRES    IN   EACH    IACRE(I)    TO   INDIVIDUAL  ACRES 
JK   =  0 

DO   30  J=l,180 

IF    (JK    . GE .    LAND)    GO   TO  4 

IF    II  ACRE ( J  )    .LE.    0)    GO   TO  30 

KL  =   JK   +  1 

JK   =   JK   +  IACRE(J) 

DO   3  I=KL,JK 

NAC   =   LAND  +1-1 

ACAGE ( NAC )    =   J   -  1 

3  CONTINUE 
30  CONTINUE 

GET   DISTRIBUTION  OF   ACRES   BY  AGE. 

CHECK   THAT  NO  ACRE    IS   OLDER   THAN   179  YEARS 

4  DO  7  K=l , LAND 

IF    ( ACAGE ( K )    .LE.    179.0)    GO   TO  6 
WRITE  (6,5) 

5  FURMAT    ( 1H1 ,47X,38HY0U  WENT   BEYOND  AGE   LIMIT   OF    179  YEARS) 
I  YEAR   =   NOYRS   -  1 


RETURN 

LM  =   AC AGE ( K )    +  1.0 

N ACRE ( L  M I    =   NACRE ( LM )    +  1 

COMPUTE   TOTAL   ACREAGE   BY    10-YEAR    AGE  CLASSES 

DO   8  1=1,18 

DO   8   J  =  l,  10 

NS  =  10  »  (I  -  11   +  J 

ISUM(I)    =    ISUM(I)    +   NACRE ( NS ) 


DO    10   M= 1 , L AND 

IF    ( ACAGE ( M )    .LT.    AGMRCH)    GO  TO  10 
I  AG   =   ACAGE ( M )    +  1.0 

IF    ( ANBDF ( I  AG )    .GE.    BFMRCH)    GO   TO  9 
GVLCU  =   GVLCU   +  ANCUV(IAG) 
GO   TO  10 
9  GVLBF   =   GVLBF   +   ANBDF (  I  AG ) 
10  CONTINUE 

COMPUTE    INITIAL   NON-ZERO  VALUES   FOR  0UTPUT2 


IVAR17, 1  ) 
-  IVAR(8,1) 


=  GVLCU  +  0. 
=  GVLBF  +  0. 
=  IVAR(5,1) 
)  =  I VARI6, 1 ] 
)  =  NONSTK 
=  PRICF(l) 
=   PRIBDI 1 ) 

GVLBF   *    (PRIBD(l)    -  BFCST) 
(GVLCU/100.0)    »    (PRICF(l)  - 
=   GSVALC   +  GSVALB 
=   VARI 13, I  )    +   VAR( 12, 1 ) 


IVARI7.1 
IVARI8, 1 
IVARI9, 1 
I VAR ( 10 , 
IVARI 11, 
VAR( 1  ,1  ) 
VAR(2,1) 
GSVALB  = 
GSVALC  = 
VAR( 13, 1 
VAR ( 14, 1 
DO    11  1=1,14 
N   =    I    +  11 
IVAR(N.l)    =  ISUM(I) 
IVAR(26,1)    =  ISUMI15 


WRITE   HEADINGS   FOR   TABLE   OF    INITIAL   DISTRIBUTION   OF   ACRES   BY  AGE 
WRITE  (6,12) 

12  FORMAT    ( 1H1,////,38X,36HINITIAL   DISTRIBUTION  OF   ACRES   BY  AGE) 
WRITE    (6,13)    (BATCH! I ), 1=1 ,3) 

13  FORMAT    ( 1H    ,45X,7HBATCH  ,3A8) 
WRITE    (6,14)  ITEST 

14  FORMAT    (1H    ,45X,4HTEST, 14) 
WRITE    (6,15)    ( GMNAM ( I ) , I = 1 , 3 ) 

15  FORMAT    ( 1H    ,45X,6HGAME  ,3A8) 
WRITE    (6,16)    (DESCRl  I  ), 1  =  1,5) 

16  FORMAT    (1H  ,45X,5A8) 
WRITE  (6,17) 

17  FORMAT    (1H    ,45X,16HYEAR  WITHIN  G AME , 3X , 1  HO , / // ) 
WRITE  (6,18) 

18  FORMAT    (1H    , 55X , 9HAGE ( YE AR ) ) 
WRITE  (6,19) 

19  FORMAT    (1H    , 4X, 1 1HAGE ( DECADE ), 7X , 1H0, 7X , 1H1 , 7X , 1H2.7X, 1H3,7X,1H4, 7 
IX, 1H5,7X,1H6,7X,1H7,7X, 1H8.7X, 1H9 , 1  OX , 5HT0T AL , //) 


WRITE  NUMBER  OF  ACRES  IN  EACH 
10-YEAR  CLASSES 


-YEAR   AGE   CLASS   AND   TOTALS  OF 


DO  21  J=l 
IK  =   J  - 


NN  =    10  *    IK   +  1 

WRITE  (6,20)  I K , NACRE ( NN ) ,  NACRE (NN+1) , NACRE (NN+2) , NACRE ( NN  +  3 ) , NACR 
1EINN+4) , NACRE (NN+5) , NACRE ( NN+ 6 ), NAC RE (NN+7), NACRE (NN+ 8), NACRE (NN+9 
2) ,  I  SUM ( J ) 

20  FORMAT    ( 1H    , I  1 1 , 5X , 10  I  8 ,  I  1 5 , / ) 

21  CONTINUE 
RETURN 
END 

SUBROUTINE  YEARS 

COMMON   ACCST, AGEO, AGMRCH, ANBDF (  1 8 1 ) , ANCU V ( 1 8 1 ) , ANNET , BATCH  (  3 )  ,  BDFC 
1(180) ,BDFO( 180) , BFCST, CFMC (180) ,CFMO( 180) , C ST AC , C ST VL , CUCST , CUT AGE 
2,DBH0,DEN0,0ESCR(5) , DL EV , FMRCHD ( 1 0 ) , GMN AM ( 3 ) , GS V AL B , G S VALC , G VL BF , 
3GVLCU, I  ACRE ( 180 ) , IALCUT, I  GAME , I  SUM ( 18) , I TE ST , I VAR ( 26 , 1 50 ) , I  YEAR, 
4K0L ( 6) , LAND , LAST, MALCUTI 10), MOLD, NACRE ( 180 ),NGAME,NK0LS,N1, NONSTK, 
5 NOYRS, PRET, PRIBDI  1 50 ) , PR  I CF (  1 50 ) , PR  I D I  V  I  10 ) , RE TRN , ROT A , RA TE , GROW , 
6SITE,SUMM(6,25, 10) , TCOS T , SHE LT , YRLOS , F I NL , CLOS S , C THN , C PLT , CFPCT , 
7VAR(14,150),VLBF,VLCU,ACAGE(1000) ,CYCL, KOUNT, DEFOR, ANUL, MIX, IPLNT, 
8RINT, THIN, BFMRCH, BFSALV, COMCU,COMBF,BFPCT, GNTR 

GVLBF   =  0.0 

GVLCU   =  0.0 

LOSS   =  0 

NPLNT   =  0 

RETHV   =  0.0 

RETTH   =  0.0 

SCLOSS  =  0.0 

SCPLT   =  0.0 

SCTHN  =  0.0 

VLBF   =  0.0 

VLCU  =  0.0 

JCYCL   =  CYCL 

IYRM  =    I  YEAR  +  1 


MAKE 


SCHEDULED   ANNUAL  PLANTING 


IF    (NONSTK    .EQ.    0)    GO  TO  1 
NPLNT   =  IPLNT 

IF    (NPLNT    .GT.    NONSTK)    NPLNT  = 

NONSTK   =  NONSTK  —NPLNT 

APLT   =  NPLNT 

SCPLT   =   APLT   *  CPLT 

IF    (NONSTK    .EQ.    0)    GO  TO  1 

KIM  =   LAND  -  NONSTK   +   1    +  LAST 
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IF    (KIM    .GT.    LAND)    KIM  =   KIM  -  LAND 
ICH   =   KIM  +   NONSTK   -  1 

C 

C  AGE    ZERO   OF   AREAS   SUBROUTINE   GIVES   CONDITIONS   AT   END   OF   YEAR  ZERO, 

C  SO    INCREASE    AGES   ONE   YEAR   TO   START  SIMULATION 

C 

1  DO   2    1=1, LAND 

ACAGE ( I )    =   AC AGE ( I )    +  1.0 

2  CONTINUE 

IF    (NONSTK    . EO .    0)    GO   TO  4 

C 

C  SUPPRESS   AGE   INCREASE  FOR  NONSTOCKED  ACRES 

C 

DO  3  [=KIH, ICH 
ACAGE ( I )    =  0.0 

3  CONTINUE 

4  IF    (DEFOR    .EO.    0.0)    GO   TO  15 

C 

C  DETERMINE   AREA   DEFORESTED  ANNUALLY 

C 

AKDX   =   LAND   -  NONSTK 

YRLOS   =    (AKDX   *   DEFOR)    +  YRLOS 

IF    (YRLOS    .LT.    1.0)    GO   TO  15 

C 

C  GENERATE   PSEUDORANDOM   NUMBER   FOR    AGE   OF    ACRE  DESTROYED 

C 

5  NDIV   =    (17.0   *    ANUL    +  3.01/128.0 

NULL   =  ANUL 

NULL    =    (17    *   NULL   +   3)    -   128    *  NDIV 
ANUL   =  NULL 

C 

C  CHECK   THAT    AGE   EXISTS   AND    IS   BETWEEN   ONE   AND   OLDEST   CURRENT  AGE 

C 

IF    (ANUL    . L  E .    0.0)    GO   TO  5 

IF    (ANUL    .GT.    ACAGE ( KOUNT ) )    GO   TO  5 
DO   6   M= 1 , LAND 
KACR  =  M 

IF    ( ACAGE ( M )    . EO.    ANUL)    GO   TO  7 

6  CONTINUE 
GO   TO  5 

C 

C  SET   LOSS   TO  REDUCE  CURRENT   ALL0WA8LE  CUT 

C 

7  LOSS   =   LOSS   +  1 
NONSTK   =   NONSTK   +  1 
YRLOS   =   YRLOS  -  1.0 

C 

C  SALVAGE   VOLUME    IF   NOT   LESS   THAN  BFSALV 

C  NO  SALVAGE   OR  CLEANUP    IF   AGE  LESS  THAN  AGEO. 

C 

IF    (NULL    .LE.    Nl)    GO   TO  9 

IF    (IYEAR    . EO.    1)    MTHN   =  FMRCHDll) 

NULL   =   NULL   +  1 

KULL   =    NULL   -  1 

IF    (KULL    .LT.    MTHN)    GO   TO  8 

SALVB    =   ANBDF ( NULL )    +  BDFC(KULL) 

IF    (SALVB    .LT.    BFSALV)    SCLOSS   =   SCLOSS   +  CLOSS 
IF    (SALVB    .LT.    BFSALV)    GO   TO  9 
VL1F   =   VLBF   +  SALVB 

RETTH   =   RETTH   +   SALVB   *    IPRIBDl IYRM)    »  BFPCT) 
GO   TO  9 

8  SALVB   =   ANBDF ( NULL ) 

IF    (SALVB   .LT.    BFSALV)    SCLOSS  =  SCLOSS   +  CLOSS 
IF    (SALVB    .LT.    BFSALV)    GO   TO  9 
VLBF   =   VLBF   +  SALVB 

RETTH   =   RETTH   +    SALVB   »    (PRIBD(IYRM)    »  BFPCT) 

C 

C  RENUMBER  ACRES   TO  PUT  ACRE   LOST   AT   END  OF   AGE  SEQUENCE 

C  WITH  AGE  ZERO 

C 

9  IF    I KACR    .NE.    KOUNT)    GO   TO  10 
LAST   =   LAST    +  1 

KOUNT  =  KOUNT  +  1 
ACAGE(LAST)  =  0.0 
GO   TO  15 

10  LUB   =    LAST   -  1 

IF    (KACR    .LT.    LAST)    GO   TO  13 
MNO   =   LAND   -  KACR 
DO    11  J=1,MN0 
JSUB  =  KACR  +  J 
ISUB  =   JSUB  -  1 
ACAGE(ISUB)    =  ACAGE(JSUB) 

11  CONTINUE 

AC AGE ( L  AND )    =   ACAGE ( 1 ) 

DO    12  K=1,LUB 

KAN  =  K  +  1 

ACAGE ( K )    =   ACAGE ( KAN ) 

12  CONTINUE 

ACAGE ( L AST )    =  0.0 
GO  TO  15 

13  DO    14   M  =  K  ACR , LUB 
MOL   =   M   +  1 

ACAGE ( M )    =    ACAGE ( MOL ) 

14  CONTINUE 

ACAGE ( LAST )    =  0.0 

IF    (YRLOS    .GE.    1.0)    GO   TO  5 

C 

C  PREPARE   SUBTOTALS   FOR  CURRENT  YEAR  AND  CHECK  THAT  NO  ACRE 

C  IS  OLDER  THAN   179  YEARS 

C 

15  DO    16  K=l,180 
NACRE (K)    =  0 

16  CONTINUE 

DO    19   K= 1 , LAND 

IF    (ACAGE(K)    .LE.    179.0)    GO   TO  18 
WRITE  (6,17) 

17  FORMAT    ( 1HI ,/////, 47X.38HY0U  WENT   BEYOND   AGE   LIMIT   OF    179  YEARS) 
IYEAR   =  NOYRS 


GO   TO  42 

18  LM   =   ACAGE ( K )    +  1.0 
NACRE ( LM )    =   NACRE ( LM )    +  1 

19  CONTINUE 

C 

C  DETERMINE   ALLOWABLE   CUT   BASED   ON   BD.FT.    STUMPAGE  PRICE 

C 

DO   20  J=l,10 
NSUB  =  J 

IF    (PRIBD(IYRM)    .LE.    PRIDIV(J))    GO   TO  21 

20  CONTINUE 

21  IALCUT   =   MALCUT ( NSUB )    -  LOSS 
CUTAGE   =   FMRCHD ( NSUB ) 

C 

C  COMPUTE   THINNINGS   FOR   ANNUAL  CUT 

C 

MX  Y  =  0 

MAC   =  CUTAGE 

DO   24    I=N1,MAC, JCYCL 

VBTH  =  0.0 

VCTH   =  0.0 

IF    (I    .GE.    MAC)    GO  TO  25 
MR   =    I    +  1 

IF    (BDFC(I)    .LT.    COMBF)    GO   TO  22 
VLBF  1   =   NACRE ( MR )    «  BDFC(I) 
VLBF   =   VLBF   +  VLBF1 
VBTH  =   VLBF 1 

RETTH   =   RETTH   +   VBTH   *    (PRIBD(IYRM)    ♦  BFPCT) 
MX Y  =   MX Y  +  1 
GO   TO  24 

22  IF    (CFMC(I)    .LT.    COMCU)    GO   TO  23 
VLCU1   =   N ACRE ( MR )    *  CFMC(I) 
VLCU   =   VLCU   +  VLCU1 

VCTH   =  VLCU1 

RETTH   =   RETTH   +   VCTH/100.0   *    (PRICF(IYRM)    »  CFPCT) 
MX Y   =   MXY  +  1 

GO   TO  24 

23  MXY  =   MXY   +  1 

SCTHN   =   NACRE  I  MR )    *   CTHN    +  SCTHN 

24  CONTINUE 

25  MTHN   =   Nl    +   MXY   *  JCYCL 

IF    I IALCUT    .LE.    0)    GO   TO  33 

C 

C  COMPUTE    VOLUME   OF   ACRES  HARVESTED 

C 

KYR   =    IYEAR   +    1    +    F I NL 

SHWD   =   SHELT   +    (SHELT   *  GROW) 

ISHWD  =  SHWD  +  0.5 

DO   32    1=1, IALCUT 

VBHV   =  0.0 

VCHV   =  0.0 

IF    (LAST    .LT.    LAND)    GO   TO  26 
LAST   =  0 

26  LAST   =   LAST   +  1 

IF    (ACAGE(LAST)    .GE.    CUTAGE)    GO   TO  27 
LAST   =   LAST  -  1 
GO   TO  33 

27  M   =   AC AGE ( L  AST ) 
K  =   M  +  1 

KOUNT   =   KOUNT   +  1 
I  SAFE   =   LAND  +  1 

IF    (KOUNT    .GE.    ISAFE)    KOUNT   =  1 

IF    ( M    .LT.    MTHN)    GO   TO  28 

VLBF2   =    ANBDF ( K )    +   BDFC(M)    -  SHELT 

IF    (VLBF2    .LT.    COMBF)    GO   TO  29 

VLDF   =   VLBF   +  VLBF2 

VBHV  =  VLBF2 

RETHV  =  RETHV  +  VBHV  *  PR  I BD (  IYRM) 
IVAR(4,KYR)  =  IVAR(4,KYR)  +  ISHWD 
VAR(3,KYR)    =   VARI3.KYR)    +    (SHWD   *   PR  I BD ( KYR ) ) 

VA  R ( 7 , K  YR )  =  VAR ( 7 , KYR )  +  SHWD  *  (BFCST  +  BFCST  «  RATE  *  FIND 
GO  TO  31 

28  VLBF2   =   ANBDF ( K )    -  SHELT 

IF    (VLBF2    .LT.    COMBF)    GO   TO  29 
VL p  c        VL^F   +  VLBF2 

vb:..  -  vLBi-^ 

ET  -1  v  =  RETHV  +  VBHV  ♦  PRIBDI  IYRM) 
VAR i 4 , KYR )  =  IVAR(4,KYR)  +  ISHWD 
;AR(3,KYR)    =   V AR ( 3 , KYR )    +■    (SHWD   *   PR  I BD ( K YR )  ) 

VAR ( 7 , KYR )  =  VAR ( 7 , KYR )  +  SHWD  »  (BFCST  +  BFCST  »  RATE  *  FIND 
GO  TO  31 

29  IF    (M    .LT.    MTHN)    GO   TO  30 
VLCU2   =   ANCUV(K)    +  CFMC(M) 

IF    (VLCU2    .LT.    COMCU)    GO   TO  32 
VLCU   =   VLCU   +  VLCU2 
VCHV  =  VLCU? 

RETHV   =   RETHV   +   VCHV/100.0   *  PRICF(IYRM) 
GO  TO  31 

30  VLCU2   =  ANCUV(K) 

IF    ( VLCU2    .LT.    COMCU)    GO   TO  32 
VLCU   =    VLCU    +  VLCU2 
VCHV  =  VLCU2 

RETHV   =   RETHV   +   VCHV/100.0   *  PRICF(IYRM) 

31  ACAGEI LAST)    =  0.0 

32  CONTINUE 

C 

C  COMPUTE  GROWING  STOCK  VOLUME.   USE   CU.FT.    IF  VOLUME    IS  LESS 

C  THAN   BFMRCH  MBF 

C 

33  DO   37    1=1, LAND 

IF    (AC.GE(I)    .LT.    AGMRCH)    GO   TO  37 
I  AG   =  ACAGEI  I  )    +  1.0 
IBG   =    I  AG  -  1 

IF    (IBG    .LT.    MTHN)    GO   TO  34 
GBL1    =   ANBDF ( I  AG  I    +  BDFC(IBG) 
IF    (GBL1    .LT.    BFMRCH)    GO   TO  35 
GVLBF   =   GVLBF   +  GBL1 
GO  TO  37 
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34  GBL1   =   ANBDFI I  AG ) 

IF    (GBL1    .LT.    BFMRCH)    GO  TO  35 
GVLBF   =   GVLBF   +  GBL1 
GO   TO  37 

35  IF    (IBG    .LT.    MTHN )    GO  TO  36 
GCL1   =   ANCUV(IAG)    +  CFMC(IBG) 
GVLCU   =   GVLCU  +  GCL1 

GO   TO  37 

36  GCL1   =   ANCUVI I  AG ) 
GVLCU   =   GVLCU   +  GCL1 

37  CONTINUE 

PREPARE   FOR  NEW   TOTALS   AND  SUBTOTALS 

DO   38  K=l,180 

38  NACRE ( K )    =  0 
DO   39  1=1,18 

39  ISUM(I)    =  0 

DO  40  K= 1 ,  L  AND 

L  M  =    AC  AGE ( K  )    +  1.0 

NACREILM)    =   NACREILM)    +  1 

40  CONTINUE 

COMPUTE   TOTAL   ACREAGE   BY    10-YEAR   AGE  CLASSES 

DO  41  1=1,18 
DO  41  J=l,10 


REMAINING  SUBROUTINES 


ISUMI I )  = 
1  CONTINUE 


COMPUTE  VOLUMES  AN6  VALUES  AT 
TO  ANUAL 


END  OF   CURRENT   YEAR   FOR  TRANSFER 


42  RETRN  =  RETTH  J 
CSTAC  =  LAND  * 
CSTVL  =  CUCST  < 
TCOST  =  CSTAC  H 
ANNET  =  RETRN  - 
GSVALB  =  GVLBF 


RETHV 

ACCST   +   SCPLT   +   SCTHN  + 
(VLCU/100.0)   +  BFCST  « 
CSTVL 
TCOST 

♦    (PRIBDI IYRM)    -  BFCST) 


SCLOSS 
VLBF 


GSVALC   =    (GVLCU/100.0)    *  (PRICF(IYRM) 


INCREASE   COSTS   ANNUALLY,    IF  DESIRED 


+    (ACCST   »  RATE) 
+    (BFCST   ♦  RATE) 
+   (CLOSS  ♦  RATE) 
(CPLT   *  RATE) 
(CTHN   «  RATE) 
+    (CUCST   »  RATE) 


ACCST  =  ACCST 
BFCST  =  BFCST 
CLOSS  =  CLOSS 
CPLT  =  CPLT  + 
CTHN  =  CTHN  + 
CUCST  =  CUCST 
RETURN 
END 

SUBROUTINE  0UTPUT1 

COMMON   ACCST, AGEO, AGMRCH, ANBDFI  1 8 1 )  ,  ANCU V ( 1 8 1 ) , ANNE T , B A TCH ( 3 ) , BDF C 
1 ( 180) ,BDFO(  180) , BFCST, CFMCI 180) ,CFMO( 180) , CSTAC , C ST VL , CUCS T , CUT  AGE 
2,DRH0,DEN0,DESCR(5) ,DLEV,FMRCHD( 10) , GMN AM ( 3), GSVALB, GSVALC, GVLBF, 
3GVLCU, I  ACRE! 180 ) , I ALCUT , I  GAME , ISUM( 18) , I TE ST , I VAR ( 26 , 1 50 )  , I  YEAR, 
4K0LI 6) , LAND, LAST, MALCUTI 10) .MOLD, NACRE ( 180) , NG 4ME , NKOL S , N 1 , NON STK , 
5N0YRS.PRET, PRIBDI 150),PRICF(150),PRIDIV(10) , RE TRN , ROT A , R A TE , GROW , 
6SITE,SUMM(6,25, 10) , TCOST , SHE LT , YRLOS , F I NL , C LOS S , CTHN , C PLT , CF PCT , 
7VAR(14,150)  , VLBF, VLCU, AC AGE  I 1000) , C YCL , KOUNT , DEFOR , ANUL , M I  X , I PLNT , 
SRI  NT, THIN, BFMRCH, BFSALV,COMCU,COMBF,BFPCT,GNTR 

WRITE   TABLE  HEADINGS 

WRITE  (6,1) 

1  FORMAT    I 1H1,////,46X,28HDISTRIBUTI0N   OF   ACRES   BY  AGE) 
WRITE    (6,2)    (BATCH (I) ,1=1,3) 

2  FORMAT    (1H    ,45X,7HBATCH  ,3A8) 
WRITE    (6,3)  ITEST 

3  FORMAT    (1H    , 45X , 4HT EST ,  I  4  ) 
WRITE    (6,4)    (GMNAMI I) , 1=1, 3) 

4  FORMAT    (1H    ,45X,6HGAME  ,3A8) 
WRITE    (6,5)    (DESCRI I ), 1=1, 5) 

5  FORMAT    (1H  ,45X,5A8) 
WRITE    (6,6)  IYEAR 

6  FORMAT    (1H    ,45X,16HYEAR   WITHIN   GAME, 16,//) 
WRITE  (6,7) 

7  FORMAT    ( 1H   , 55X , 9HAGE ( YE AR ) ) 
WRITE  (6,8) 

8  FORMAT    (1H  ,4X,11HAGE(DECADE),7X,1H0,7X,1H1,7X,1H2,7X,1H3,7X,1H4,7 
1X,1H5,7X,1H6,7X,1H7,7X,1H8,7X,1H9,10X,5HT0TAL,//) 

WRITE   NUMBER  OF   ACRES    IN   EACH   1-YEAR   AGE   CLASS   AND   TOTALS  OF 
10-YEAR  CLASSES 

DO   10  J=l,18 

K  =  J  -  1 

NN  =    10   *  K   +  1 

WRITE  (6,9)  K , NACRE ( NN ) , NACRE ( NN+ 1 ), NACRE < NN+2 ), NACRE ( NN+3 ), NACRE < 
1NN+4) , NACRE ( NN+5 ) , NACRE ( NN+6 ),NACRE(NN+7),NACRE(NN+8),NACRE(NN+9), 

?  I  SUM ( J ) 

9  FORMAT    (1H    ,111, 5X, 1018, 115,/) 
10  CONTINUE 

RETURN 
END 

SUBROUTINE  ANUAL 

CUfMON   ACCST, AGEO, AGMRCH, ANBDFI  18  1) , ANCUVI  1 8  1  ) , ANNET , BATCH ( 3 ) , 8DFC 
1( 180) ,BDFO( 180) , BFCST, CFMCI 180) ,CFMO( 180 ), CSTAC , CSTVL , CUCST , CUT AGE 
2,DBHn,DEN0,DESCR(5),DLEV,FMRCHD(10) , GMNAMI 3 ), GSVALB, GSVALC, GVLBF, 
3GVLCU, I  ACRE  I  180) , I ALCUT , I  GAME, I  SUM  I  18) , ITEST, IVARI 26,  150) , IYEAR, 
4K0L (6) , LAND, LAST, MALCUTI 10) , MOLD, NACRE! 180) ,NGAME,NK0LS,N1 »NONSTK, 
5N0YRS,PRET,PRIBD( 150) , PR  I CF ( 150 ) , PR  I D I V ( 10 ) , RETRN , ROT A , R ATE , GROW , 
6S I TE , SUMM (6, 25, 10), TCOST, SHELT, YRLOS, FINL, CLOSS, CTHN, CPLT, CFPCT, 
7 VAR I  14,150)  , VLBF, VLCU, AC AGE (1000) , C YCL , KOUNT ,  LIE  FOR , ANUL , M I  X , I PLNT, 
HRI NT, THIN, BFMRCH, BFS AL V , COMCU , COMBF , BFPCT.GNTR 

CONVERT   VOLUME   AND   AREA   VALUES   TO   SUBSCRIPTED   VALUES   FOR  USE  BY 


K  =  IYEAR 
J  =  IYEAR  + 
IVAR(l.J)  = 
IVARI2.J)  = 
IVARI3.J)  = 
IBFT  =  VLBF 
IVARI4.J)  = 
IVAR(5,J)  = 
IVAR(6,J)  = 
IVARI7.J)  = 
IVAR(8,J)  = 
IVARI9.J)  = 
IVARIIO.J)  = 
IVARIll.J)  = 
DO  1  1=1,14 
N  =    I    +  11 


I ALCUT 

CUTAGE 

VLCU  +  0.5 

+  0.5 

I VAR I  4 , J ) 

IVARI5.K) 

IVARI6.K) 

GVLCU   +  0. 

GVLBF   +  0. 

I VAR ( 5 , J ) 

IVARI6.J) 
■  NONSTK 


IBFT 
IVARI 3, 
IVARI4, 


I VAR I  7 , J ) 
•  IVARI8.J) 


1    IVAR(N.J)    =  ISUMII) 
IVARI26.J)    =  ISUMI15) 


ENTER   MONEY   VALUES    IN   ARRAYS   FOR   REMAINING  SUBROUTINES 


VAR(1,J) 
VAR I  2, J) 
VAR(3,J) 
VAR ( 4 , J ) 
VAR(5,J) 
VAR I  6 , J ) 
VAR I  7 , J ) 
VAR I  8 , J ) 
VAR I  9 , J ) 
VARI 10, J ) 
VARIll, 
VAKI 12, 
VARI 13, 
VARU4, 
RETURN 
END 

SUBROUTINE  0UTPUT2 

COMMON  ACCST, AGEO, AGMRCH, ANBDFI 181), ANCUV (181), ANNET , BATCH  I  3 ) , BDFC 
11180) ,BDFO( 180) , BFCST, CFMCI 180 ) ,CFMO( 1 80 ) , C STAC , C STVL , CUC ST , CUT AGE 
2 , DBHO ,DEN0,DESCR(5),DLEV,FMRCHD( 10) , GMNAMI 3) , GSVALB , GS VALC , GVLBF , 
3GVLCU, I  ACRE  I  180) , I ALCUT , I  GAME , I  SUM ( 18 ), I  TEST , I VAR I  26 ,  150), IYEAR, 
4K0L (6 ) , LAND, LAST, MALCUTI 10), MOLD, NACRE (180) , NGAME , NKOL S , Nl , NONSTK , 
5N0YRS,PRET,PRIBD( 150) , PR  I C F I  150 ) , PR  I D I V  I  1 0 ) , RETRN , ROT A , R ATE , GROW , 
6SI TE, SUMM I  6, 25, 10 ), TCOST, SHELT, YRLOS, FINL, CLOSS, CTHN, CPLT, CFPCT, 
7VAR1 14, 150) ,VLBF, VLCU, AC  AGE  I  1000)  , C YCL , KOUNT , DE F OR , ANUL , M I  X , I PLNT , 
8R I  NT, THIN, BFMRCH, BFS ALV, COMCU, COMBF, BFPCT.GNTR 

PRINT   FIRST  PAGE 


PRICFI J) 
PRIBDI J ) 
VAR (3, J)  ■ 
VARI4.K)  • 
CSTAC 
V  AR  (  6  ,  K  )  • 
VAR I  7 , J )  • 
VARI8.K)  - 
TCOST 
=  VARI10.K 
=  ANNET 
=  VARI12.K 
=  GSVALC  + 
=   VARI 12, J 


VARI5, J) 

CSTVL 

VAR(7,J) 

+    V AR ( 9 , J ) 

+   VAR ( 11, J) 
GSVALB 
+  VARI 13, J  ) 


NOYRS  +  1 


1  DO  15 
LINE  = 

2  IF    I M 


.LT.   40)    GO   TO  11 


3A8) 


WRITE    (6,3)    (BATCH! I ), 1=1,31 

3  FORMAT    ( 1H1,//,46X,7HBATCH 
WRITE    (6,4)  ITEST 

4  FORMAT    ( 1H   , 45X , 4HTE ST , I  4 ) 
WRITE    (6,5)    (GMNAMI I ), 1=1, 3) 

5  FORMAT    (1H    ,45X,6HGAME  ,3A8) 
WRITE    (6,6)    (DESCRI I ), 1=1,5) 

6  FORMAT    (1H  ,45X,5A8) 
WRITE  (6,7) 

7  FORMAT    (1H  ,/) 
WRITE  (6,8) 

8  FORMAT    ( 1H    , 1 2X , 9HALL0WA BL E , 5X , 7HCUTT I NG , 8X , 10H ACTUAL  CUT,10X,9HCU 
1MUL   CUT,10X,9HGRSTK   VOL , 1 2 X , 9HT0T AL  VOL) 

WRITE  (6,9) 

9  FORMAT    (1H    , 2X , 4H YE AR , 9X , 3HCUT , 10X , 3HAGE , 7X , 6HCU . FT . , 5 X , 3HMBF , 6X , 6 
1HCU.FT.,5X,3HMBF,6X,6HCU.FT., 5X , 3HMBF , 6X , 6HCU. F T . , 5  X , 3HMBF ) 

WRITE  (6,10) 

10  FORMAT    (1H    , 1 5X , 3H ( 1 ) , 1  OX , 3H ( 2 ) , 8X , 3H ( 3 ) , 7X , 3H ( 4 ) , 7X, 3H ( 5 ) , 7X , 3H ( 6 
1 ) ,7X, 3H( 7) ,7X,3H(8) ,7X,3H(9) ,6X,4H( 10) ,//) 

11  WRITE    (6,12)    LINE,  ( IVARI I, J)  ,  1  =  1,10) 

12  FORMAT    I 1H  ,16,112,113,112,19,3(111,19)) 


IF    I  J    .LE.    1 )    GO   TO  13 

M  =   M  +  1 

IF    ILL    .LT.    10)    GO   TO  14 


13  WRITE  (6,7) 
LL   =  0 

14  LL   =   LL   +  1 

15  CONTINUE 

PRINT   SECOND  PAGE 

M  =  40 

DO  23  J=1,N 

LINE  =  J  -  1 

IF    (M   .LT.    40)    GO  TO  19 

M  =  0 

WRITE    (6,3)    (BATCH! I ), 1=1,3) 

WRITE    (6,4)    I  TEST 

WRITE    (6,5)    (GMNAMI I ) ,1=1,3) 

WRITE    (6,6)    (DESCRI I ), 1=1,5) 

WRITE  16,7) 

WRITE  (6,16) 

16  FORMAT    (1H    , 1 1 X , 3HN0N , 46X , 1 1H AGE  CLASSES) 
WRITE  (6,17) 

17  FORMAT    (1H    ,2X,4HYEAR,5X,3HSTK, 114H        0-9      10-19      20-29     30-39  40 
1-49      50-59     60-69      70-79      80-89     90-99      100-109      110-119  120-129 
2    130-139  140-179) 

WRITE  (6,18) 

18  FORMAT    (1H    , 1  OX , 1 1 7H I  1 1 )         (12)-      (13)         (14)        (15)         (16)  (17) 
1      (18)         (19)        (20)        (21)  (22)  (23)  (24)  (25) 


21 


2  (26),//) 

19  WRITE    (6,20)    L I NE , ( I VAR ( I , J ) , I = 1 1 , 26 ) 

20  FORMAT    (1H  ,16,217,16,817,18,419) 
IF    (J    .LE.    1)    GO   TO  21 

M   =   M  +  1 

IF    (LL   .LT.    10)    GO  TO  22 

21  WRITE  (6,7) 

LL  =  0 

22  LL   =    LL   +  1 

23  CONTINUE 

PRINT    THIRD  PAGE 

H  =  40 

00   31  J=1,N 

LINE  =   J  -  1 

IF    (M   .LT.    40)    GO   TO  27 

M  =  0 


WRITE    (6,3)    (BATCH! J ) ,1=1,3) 

WRITE    (6,4)  ITEST 

WRITE    (6,5)    (GMNAMI I ) , 1=1 ,3) 

WRITE    (6,6)  (DESCR(I),I=1,5) 

WRITE  (6,7) 

WRITE  (6,24) 

24  FORMAT    (1H    , 1RX , 14HSTUMPAGE   PRICE, 9X, 15HSTUMPAGE    I NCOME , 1 3X , 10HARE 
1A  COSTS, 15X,12H VOLUME  COSTS) 

WRITE  (6,25) 

25  FOkMAT    (1H   , 2X , 4HYE AR , 9X , 1 0H1 00   CU. F T ., 5X , 3HMBF , 6X , 6H ANNUAL , 5X , 9HC 
1UMULATED.6X , 6H ANNU AL , 5 X , 9HCUMUL AT  ED , 6X , 6HANNU AL , 5X , 9HCUMUL ATED ) 

WRITE  (6,26) 

26  FORMAT    ( 1H    , 1 8X , 4H ( 27 ) , 8 X , 4H I  28 ) ,6X,4H(29) ,9X,4H( 30) ,9X, 4H ( 31 ) , 9X , 
14H( 32) ,9X,4H(33 ) ,9X,4H(34) ,/  /  ) 

27  WRITE    (6,28)    LINE , ( VAR ( I , J ) , I  =  1 , 8  ) 

20  FORMAT    (1H    , 1 6 , F 1 6 . 2 , F 1 1 . 2 , F 1 2 . 0 , F 1 1 . 0 , 2 ( F 1 5 .0 , F 1 1 . 0 ) ) 
IF    (J    .LE.    1)   GO  TO  29 

IF    (LL    .LT.    10)    GO   TO  30 
29   WRITE  (6,7) 


PRINT   FOURTH  PAGE 

M  =  40 

DO    39  J=1,N 

LINE  =  J  -  1 

IF    (M    .LT.    40)    GO   TO  35 

M  =  0 

WRITE    (6,3)  (BATCH(I),I=1,3) 

WRITE   (6,4)  ITEST 

WRITE    (6,5)  (GMNAM(I),I=1,3) 

WRITE    (6,6)  (DESCR(I),I=1,5) 

WRITE  (6,7) 

WRITE  (6,32) 

32  FORMAT    ( 1H   ,  18X , 10HT0TAL  COST , 17X , 10HNET   I  NCOME , 1 3X , 1 3HCURR ENT  VAL 
1UE,9X,5HT0TAL( 

WRITE  (6,33) 

33  FORMAT    ( 1H    , 2 X , 4H YE AR , 8X , 6HANNU AL , 5X , 9HCUMUL AT ED , 7X , 6H ANNUAL , 5X , 9H 
1CUMULATED,7X, 13HGR OWING  STOCK, 7X,9HNET  WORTH) 

WRITE  (6,34) 

34  FORMAT    ( 1H   , 1 5X , 4H ( 35 ) , 9X , 4H ( 36 ) , 10X , 4H ( 37 ) , 9X , 4H ( 3 8 ) , 1 3X , 4H ( 39 ) , 1 
14X,4H(40) ,//) 

35  WRITE    (6,36)    LINE,(VAR(  I, J), 1  =  9, 14) 

36  FORMAT    (1H   , 1 6 , F 14 . 0 , F 1 2 . 0 , F 1 5 . 0 , F 1 2 . 0 , 2F 1 8 . 0 ) 
IF    (J    .LE.    1)   GO  TO  37 

M   =    M   +  1 

IF    (LL    .LT.    10)    GO   TO  38 

37  WRITE  (6,7) 
LL   =  0 

38  LL   =   LL  +  1 

39  CONTINUE 
RETURN 
END 

SUBROUTINE  WORTH 

COMMON   ACCST,AGEO,AGMRCH,ANBDF( 181) , ANCUVI 1 8 1 ) , ANNET , B ATCH ( 3 ) , BDFC 
1  (  180) ,BDFO(  180) ,BFCST,CFMC( 180 ) ,CFMO( 1 8 0 ) , C ST AC , C ST VL , C UC ST , CUT AGE 
?,DBHO,DENO,DESCR( 5) , OLE  V , F MRCHD ( 1 0 ) , GMN AM ( 3 ) , GS V AL B , GS V ALC , GVL BF , 
3GVLCU , I  ACRE ( 180) , I AL CUT , I G AME , I  SUM! 18) , I TE ST , I VAR ( 26 , 1 50 ) , I  YEAR, 
4K0LI6) , LAND, LAST, MALCUT( 10) .MOLD, NACRE ( 1 80 ) , NG AME , NKOL S , N 1 , NONSTK , 
5N0YRS,PRET,PRIBD( 150) ,PRICF( 150) ,PRI0IV( 10) ,RETRN, ROTA, RATE, GROW, 
6SITE,SUMM(6,25, 1 0 ) , TCOS T , SHELT , YRLOS , F I NL , CLOSS , CTHN , CPLT ,  CF PCT  , 
7VAR(  14, 150) , VLBF.VLCU, AC  AGE ( 1000) , C YCL , KOUNT , DEFOR , ANUL , M I  X , IPLMT, 
8RINT,THIN,BFMRCH,BFSALV,COMCU,C0MBF,BFPCT,GNTR 

DIMENSION  CRATE (20), DISCI 20) ,01  SGI  20) ,01  SI (20) , PREV( 20) , PWTH(20),R 
1ATI0I20, 150) 

DO   1  1=1,20 
•      CRATE ( I )    =  0.0 


DISCI  I  ) 
DISG( I  ) 
DISI  (  I  ) 
PREVI  I  ) 
PWTHI  I  ) 

DO  2  t=: 


=  0.0 
=  0.0 
=  0.0 
=  0.0 
=  0.0 
,20 


DO  2  J=l,150 
RATIO! I , J)    =  0.0 


PRF  PARE   ARRAY   OF    ALTERNATIVE  RATES 


CRATE! 1 )  =  0.010 
DO   3  1=1,19 


3   CRATF(K)    =   CRATE! I )    +  0.005 

PREPARE    INTEREST   TABLE  FOR  PERIOD  NOYRS 


DO  4  K=l, NOYRS 
YRS   =  K 

PROD  =   A LOG  I FACTR )   •  YRS 
RATIOIJ.K)    =  EXPIPROD) 

4  CONTINUE 

DISCOUNT  GROWING  STOCK  VALUE   AT  NOYRS 

DO   5  L=l,20 
KL    =   NOYRS   +  I 

DISG(L)   =  VAR! 13, KLI/RATIOIL, NOYRS) 

5  CONTINUE 

DISCOUNT   ANNUAL  COSTS  AND  RETURNS 


DO 


1,20 


N) 


PRESC  =  0.0 
PRESI  =  0.0 
SPRSC  =  0.0 
SPRSI  =  0.0 
DO  6  N=l, NOYRS 
I   =  N  +  1 

PRESC   =   VAR(9, I )/RATIO(M,N 
PRE  S I    =   VAR(3, I  1/RATIOI 
SPRSC  =  SPRSC  +  PRESC 
SPRSI   =  SPRSI   +  PRE  S I 
CONTINUE 
DISKM)    =  SPRSI 
DISC(M)    =  SPRSC 
CONTINUE 


OBTAIN   PRESENT   WORTH   AT    EACH  RATE 


DO  4  J=l,20 
FACTR   =  1.0 


DO   8  IJ=1,20 

PREVI  IJ)   =   DISI ( IJ)    +  DISGI  IJ ) 

PWTH(IJ)    =   PREVIIJ)    -   DISC(IJ)    -  VAR(13,1) 

CRATE(IJ)   =   CRATE(IJ)    •  100.0 

8  CONTINUE 
WRITE  (6,9) 

9  FORMAT    I 1H1,////,53X,29HPRESENT   WORTH   AND  RATE  EARNED) 
WRITE    (6,10)    (BATCH!  I), 1  =  1, 3) 

10  FORMAT    (1H    ,52X,7HBATCH      , 3A8 ) 
WRITE    (6,11)  ITEST 

11  FORMAT    (1H   ,52X,4HTEST, 14) 
WRITE    (6,12)    (GMNAMI I  ), 1  =  1,3) 

12  FORMAT    (1H   ,52X,6HGAME  ,3A8) 
WRITE    (6,13)    (DESCRI I ), 1=1,5) 

13  FORMAT   ( 1H  ,52X,5A8) 
WRITE    (6,14)  NOYRS 

14  FORMAT    ( 1H    , 52X , 1 5HYE AR S    IN   PERIOD, 15,//) 
WRITE    (6,15)  VAR(13,1) 

15  FORMAT    (1H   , 1 1 X , 33HV ALUE  OF   INITIAL   GROWING  STOCK  —  J , F 1 0. 2 , // ) 
WRITE  (6,16) 

16  FORMAT    (1H    ,  57X ,  38HVALUES  DISCOUNTED  TO  PRESENT  (DOLLARS),/) 
WRITE  (6,17) 

17  FORMAT   (1H   , 1 1 X , 8HC0MP0UND , 1 5 X , 6HFUTURE , 34X , 5HST0CK , 36X , 3HNET ) 
WRITE  (6,18) 

18  FORMAT    ! 1H   , 13X , 4HR ATE , 1 6X , 7HGR0W I NG , 1 5 X , 3H ALL , 1 7 X , 4HPLUS , 1 6X ,  3HAL 
1L,15X,7HPRESENT) 

WRITE  (6,19) 

19  FORMAT    (1H   , 1 1 X , 9H ( PERC ENT ), 1 4X , 5HST0CK , 14X , 7H  I NCOMES , 1 3X , 7H INCOME 
1S,14X,5HC0STS, 15X, 5HW0RTH, / ) 

DO  21  1=1,20 

WRITE    (6, 20) CRATE! I ) , D I SG (  I ) , D I S I ( I ) , PRE V ( I ) , D I SC ( I ) , P WTH (  I  ) 

20  FORMAT    (1H   ,  1 2X  ,  F5 . 1 ,  1 2X , 5 ( F  10 . 2 , 1 OX ) , / ) 

21  CONTINUE 
RETURN 
END 

SUBROUTINE  SUMRY 

COMMON  ACCST.AGEO, AGMRCH, ANBDFI 181 ) , ANCUVI 181) , ANNET , BATCH  I  3 ),  BDFC 
II 180) ,BDF0(180) ,BFCST,CFMC( 180) ,CFMO( 180)  , C ST AC , C ST VL , CUCST , CUT AGE 
2,DBH0,DEN0,DESCR(5),DLEV,FMRCHD(  1 0  I , GMN AM ( 3  )  , G S VALB , GS V ALC , G VLBF  , 
3GVLCU, IACRE! 180) , I ALCUT, I GAME, I  SUM ( 18) , I TE ST , I  V AR ( 26 , 1 50 ) , I  YEAR, 
4K0LI6) ,LAND,LAST,MALCUTI 10) , MOLD, NACRE!  1 80  )  , NG AME , NKOL S , N 1 , NONSTK , 
5N0YRS,PRET,PRIBD( 1 50 ) , PR  I CF ( 1 50 ) , PR  I D I V ( 10 ) , RE TRN  ,  ROT A , RATE , GROW , 
6SITE,SUMM(6,25, 10) , TCOST , SHELT , YRLOS , F I N L , CLOS S , CTHN , CPLT , CFPCT , 
7VAR( 14, 150) ,VLBF, VLCU, AC AGE ( 1 000 ) , C YCL , KOUNT , DE FOR , ANUL , M I  X , I PLNT , 
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF,BFPCT,GNTR 

CONVERT   IVARII.J)    AND  VAR(I,J)    TO  SUMMfl.J.K)    AT  END  OF  EACH  GAME 

LIM  =   10  +  NOYRS/10 
DO  4  I=1,NK0LS 
DO  4  J=1,LIM 
K  =   KOL I  I ) 

IF    (J   .GT.    10)    GO  TO  1 
JJ  =   J   +  1 
GO  TO  2 

1  JJ   =   10  • 

2  IF    (K  .GT 
SUMMI I , J, IGAME) 
GO  TO  4 

3  K  =   K  -  26 
SUMMI I , J, IGAME)    =  VAR(K,JJ) 

4  CONTINUE 

WRITE   SUMMARY  TABLES   IF   ALL  GAMES  FINISHED 
IF    (IGAME   .LT.   NG AME )    GO  TO  17 

WRITE   PAGE  HEADINGS.    SEPARATE  PAGE  FOR   EACH  COLUMN  OF  0UTPUT2 
IDENTIFIED   IN  INPUT1 

DO   16  I=1,NK0LS 
WRITE  (6,5) 

5  FORMAT    (1HI,////,46X,26HC0MPARIS0N  OF  ALTERNATIVES) 
WRITE    (6,6)    ( BATCH! I ), 1=1, 3) 

6  FORMAT    (1H    ,45X,7HBATCH      , 3A8 ) 


J  -  10)  +  1 
26)    GO  TO  3 

IVARIK, JJ) 


-  22  - 


ITEST 
,45X,4HTEST, 

1DESCRI I ) , 1 
,45X,5A8> 


14) 


13,///) 


WRITE  (6,7) 

7  FORMAT  (1H 
WRITE  (6,8) 

8  FORMAT    (1H  , 
K  =  KOL(I) 
WRITE    (6,9)  K 

9  FORMAT    ( 1H    , 45X , 8HCOLUMN 
WRITE  (6,10) 

10  FORMAT  (1H  ,  5X  ,  4HYEAR , 6X , 6HGAME  1,6X,6HGAME  2,6X,6HGAME  3,6X,6HG 
IE  4,6X,6HGAME  5,6X,6HGAME  6,6X,6HGAME  7 , 6X , 6HGAME  8,6X,6HGAME  9, 
2.7MGAME  10,//) 


WRITE  SUMM(I,J,K)  FOR  EACH  OF  FIRST  10  YEARS  AND  FOR  END  OF 
EACH  DECADE 


DO  16 
IF    I J 
JJ   =  J 
GO   TO  12 

11  J J   =    10   »    ( J   -  10) 

12  IF    ( H    . L T .    5)    GO   TO  14 
WRITE  (6,13) 

13  FORMAT    (1H  ,//) 
M  =  0 

14  WRITE    (6,15)    JJ, ( SUMM( I , J,L) , 

15  FORMAT    (1H    ,  I  9 , F 1 1 . 0 , 9F 1 2 . 0 ) 

16  M   =  M+l 

17  RETURN 
END 


APPENDIX  2 
Output  of  Test  Problem 


PAGE  TYPE  I 


YIELDS   PER  ACRE  OF   MANAGED,   EVEN-AGED  STANDS  OF   PONDE ROSA  PINE   IN  THE  BLACK  HILLS 


SITE 

INDEX  60, 

20-YEAR 

CUTTING  CYCLE,  DENSITY 

LEVEL 

100 

ENTIRE 

STAND  BEFORE  AND 

AFTER  THINNING 

PERIODIC  CUT 

AND  MORT  AL I T> 

STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
MBF 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

30 
30 

1000 
496 

110 

73 

4.5 
5.2 

20 
20 

818 
534 

201 
201 

0. 
0. 

504 

37 

284 

0 

40 

496 

104 

6.2 

28 

1146 

684 

0. 

50 
50 

496 
286 

133 
86 

7.0 
7.4 

35 
35 

1866 
1209 

1381 
954 

0. 
0. 

210 

47 

657 

427 

60 

286 

108 

8.3 

41 

1806 

1556 

1.790 

70 
70 

286 
200 

129 
99 

9.1 
9.5 

47 
47 

2501 
1930 

2241 
1756 

3.880 
3.510 

86 

30 

571 

485 

80 

200 

116 

10.3 

52 

2515 

2341 

5.990 

90 
90 

200 
141 

132 
100 

11.0 
11.4 

57 
57 

3144 
2407 

2940 
2255 

8.990 
7.560 

59 

32 

737 

685 

100 

141 

115 

12.2 

60 

2902 

2730 

10.590 

110 

no 

141 
102 

130 
100 

13.0 
13.4 

63 
63 

3457 
2701 

3263 
2555 

13.650 
11.020 

39 

30 

756 

708 

120 

102 

114 

14.3 

66 

3216 

3055 

13.960 

130 
130 

102 
76 

127 
100 

15.1 
15.5 

69 
69 

3741 
2998 

3565 
2863 

16.980 
13.910 

26 

27 

743 

702 

140 

76 

112 

16.4 

71 

3445 

3297 

16.740 

150 

76 

123 

17.2 

73 

3887 

3723 

19.510 

SAWTIMBER 
VOLUME 
MBF 


PAGE  TYPE  2 


GROWING  STOCK  OF  MANAGED  BLACK  HILLS  PONDEROS A  PINE 
SITE   INDEX  60,     20-YEAR  CUTTING  CYCLE 
DENSITY  LEVEL-   120  AND  100 


VOLUMES  PRESENT  PER  ACRE  AT  ENO  OF  EACH  YEAR 
MERCHANTABLE  CUBIC  FEET 


DECADE 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3 

201.0 

249.3 

297.6 

345.9 

394.2 

442.5 

490.8 

539.1 

587.4 

635.7 

4 

684.0 

753.7 

823.4 

893.1 

962.8 

1032.5 

1102.2 

1171.9 

1241.6 

1311.3 

5 

954.0 

1014.2 

1074.4 

1134.6 

1194.8 

1255.0 

1315.2 

1375.4 

1435.6 

1495.8 

6 

1556.0 

1624.5 

1693.0 

1761.5 

1830.0 

1898.5 

1967.0 

2035.5 

2104.0 

2172.5 

7 

1756.0 

1814.5 

1873.0 

1931.5 

1990.0 

2048.5 

2107.0 

2165.5 

2224.0 

2282.5 

8 

2341.0 

2400.9 

2460.8 

2520.7 

2580.6 

2640.5 

2700.4 

2760.3 

2820.2 

2880.1 

9 

2255.0 

2302.5 

2350.0 

2397.5 

2445.0 

2492.5 

2540.0 

2587.5 

2635.0 

2682.5 

10 

2730.0 

2783.3 

2836.6 

2889.9 

2943.2 

2996.5 

3049.8 

3103.1 

3156.4 

3209.7 

11 

2555.0 

2605.0 

2655.0 

2705.0 

2755.0 

2805.0 

2855.0 

2905.0 

2955.0 

3005.0 

12 

3055.0 

3106.0 

3157.0 

3208.0 

3259.0 

3310.0 

3361.0 

3412.0 

3463.0 

3514.0 

13 

2863.0 

2906.4 

2949.8 

2993.2 

3036.6 

3080.0 

3123.4 

3166.8 

3210.2 

3253.6 

14 

3297.0 

3339.6 

3382.2 

3424.8 

3467.4 

3510.0 

3552.6 

3595.2 

3637.8 

3680.4 

15 

3723.0 

-  23  - 


PAGE  TYPE  2  (continued) 


THOUSANDS   OP   BOARD  FEET 


0. 
0. 
0. 
0. 

1.790 
3.510 
5.990 
7.560 
10.590 
11.020 
13.960 
13.910 
16.740 
19.510 


0. 
0. 

.179 
1.999 
3.  758 
6.290 
7.863 
10.896 
11.314 
14.262 
14.193 
17.017 


.358 
2.208 
4.006 
6.590 
8.  166 
11.202 
11.608 
14.564 
14.476 
17.294 


0. 
0. 

c. 

.537 
2.417 
4.254 
6.890 
8.469 
11.508 
11.902 
14.866 
14.759 
17.571 


0. 
0. 
0. 

.716 
2.626 
4.502 
7.190 
8.  772 
11.814 
12.196 
15.168 
15.042 
17.848 


0. 
0. 
0. 
0. 
0. 

.895 
2.835 
4.750 
7.490 
9.075 
12.  120 
12.490 
15.470 
15.325 
18.125 


0. 
0. 
0. 
0. 
0. 

1.074 
3.044 
4.998 
7.790 
9.378 
12.426 
12.784 
15.772 
15.608 
18.402 


1.253 
3.253 
5.246 
8.090 
9.681 
12.732 
13.078 
16.074 
15.891 
18.679 


0. 
0. 

0. 
0. 

0. 

1.432 
3.462 
5.  494 
8.390 
9.984 
13.038 
13.372 
16.376 
16. 174 
18.956 


0. 
0. 
0. 

1.611 
3.671 
5.742 
8.690 
10.287 
13. 344 
13.666 
16.678 
16.457 
19.233 


PAGE  TYPE  3  alternatives  for  this  game 

BATCH      TEST  PROBLEM 
TEST  1 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGED i THINNED   AGE  30 


NUMBER   OF   YEARS   PER   GAME  30 


CRITICAL   PRICES  99.00  0.  0.  0.  0.  0.  0.  0.  0.  0. 

ALLOWABLE  CUT  7000000000 
MINIMUM   CUTTING   AGE  130  000000000 


ACRES   IN  WORKING  CIRCLE  915 

MINIMUM   VALUES   FOR    INCLUSION    IN  TOTALS 

AGE,    FOR   GROWING   STOCK  40 

M   BD.    FT. ,   FOR  GROWING  STOCK  1.5 

CU.    FT. ,    FOR   COMMERCIAL   CUT  400 

M   BD.    FT.,    FOR   COMMERCIAL   CUT  3.0 

M    BD.    FT.,    FOR   SALVAGE  1.5 


COSTS    IN   FIRST   YEAR   OF  GAME 

PER  ACRE    (ANNUAL)  .20 

PER   100  CU.    FT.    HARVESTED  .05 

PER   M   BD.    FT.  1.56 

THIN   ONE    ACRE  25.00 

PLANT    ONE    ACRE  30.00 

CLEANUP  OF   ONE   ACRE  25.00 

RATE   OF    INCREASE    IN   COSTS  .01 


ACRES    PLANTED  ANNUALLY 
PERCENT   OF   ACRES   LOST  ANNUALLY 
M   BD.    FT.    IN  SHELTERWOOD 
REGENERATION  PERIOD 


1 

.040 
4.0 
10.0 


RELATIVE    VALUE   OF    INTERMEDIATE  CUTS 

STUMPAGE  PRICE,  CU.  FT.  1.00 
STUMPAGE   PRICE,    BD.    FT.  .85 


PSEUDORANDOM   NUMBER  GENERATOR 


PAGE  TYPE  4 


INITIAL   DISTRIBUTION   OF    ACRES   BY  AGE 
BATCH      TEST  PROBLEM 
TEST  1 

GAME      EOUAL    AREAS   CUT  ANNUALLY 
MANAGED, THINNED  AGE  30 
YEAR   WITHIN   GAME  0 


AGE ( DECADE ) 
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PAGE  TYPE  4  (continued) 


DISTRIBUTION   OF    ACRES   BY  AGE 
BATCH      TEST  PROBLEM 
TEST  1 

GAME      EQUAL    AREAS   CUT  ANNUALLY 
MANAGED, THINNED   AGE  30 
YEAR   WITHIN   GAME  30 


AGE(DECADE) 


PAGE  TYPE   5  BATCH      TEST  PROBLEM 

TEST  I 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGED, THINNED   AGE  30 


ALLOWABLE 
CUT 
(  1  ) 


ACTUAL  CUT 
CU.FT.  MBF 
(3)  14) 


CUMUL  CUT 
CU.FT.  MBF 
(5)  (6) 


GRSTK  VOL 
CU.FT.  MBF 
(7)  18) 


CU.FT. 
19) 


MBF 
(10) 


130 
130 
130 
130 
130 
130 
130 
130 
130 
130 


16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
14719 
16135 


16135 
32270 
48405 
64540 
80675 
96810 
1 12945 
129080 
143799 
159934 


182 
272 
363 
454 
545 
636 
727 
805 
897 


145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 


4439 
4439 
4444 
4444 
4443 
4447 
4447 
4446 
4468 
4467 


161243 
177378 
193513 

209648 
225783 
241918 
258053 
274188 
288907 
305042 


4530 
4621 
4716 
4807 
4897 
4992 
5083 
5173 
5273 
5364 


130 
130 
130 
130 
130 
130 
130 
130 
130 
130 


16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
15450 


114 

127 
122 
114 
127 
122 
122 
128 
123 
119 


176069 
192204 
208339 
224474 
240609 
256744 
272879 
289014 
305149 
320599 


1011 
1138 
1260 
1374 
1501 
1623 
1745 
1873 
1996 
2115 


145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
144424 
144354 


4483 
4483 
4482 
4499 
4498 
4497 
4506 
4505 
4504 
4517 


321 177 
337312 
353447 
369582 
385717 
401852 
417987 
434122 
449573 
464953 


5494 
5621 
5742 
5873 
5999 
6120 
6251 
6378 
6500 
6632 


130 
130 
130 
130 
130 
130 
130 
130 
130 
130 


16135 
16135 
15427 
16135 
16135 
16135 
16135 
16135 
15708 
16135 


123 
128 
128 
123 
117 
128 
123 
118 
129 
124 


336734 
352869 
368296 
384431 
400566 
416701 
432836 
448971 
464679 
480814 


2238 
2366 
2494 
2617 
2734 
2862 
2985 
3103 
3232 
3356 


144284 
144215 
144145 
144075 
144006 
143936 
143866 
143797 
144154 
144094 


4516 
4519 
4520 
4518 
4531 
4530 
4528 
4540 
4539 
4537 


481018 
497084 
512441 
528506 
544572 
560637 
576702 
592768 
608833 
624908 


6754 
6885 
7014 
7135 
7265 
7392 
7513 
7643 
7771 
7893 
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PAGE  TYPE  5  (continued) 


BATCH      TEST  PROBLEM 
TEST  1 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGED, THINNED   AGE  30 


NON 

STK  0-9  10-19 
(11)        (12)  (13) 


20-29 
(  14) 


50-59 
(17) 


AGE  CLASSES 
60-69     70-79  80-89 
(18)        (19)  (20) 


90-99 
(21) 


100-109 
(  22) 


110-119     120-129     130-139  140-179 
(23)  (24)  (25)  (26) 


BATCH     TEST  PROBLEM 
TEST  1 

GAME  EQUAL  AREAS  CUT  ANNUALLY 
MANAGED, THINNED  AGE  30 


STUMPAGE  PRICE 
100  CU.FT.  MBF 
(27)  (281 


STUMPAGE  INCOME 
ANNUAL  CUMULATED 
(29)  (30) 


AREA  COSTS 


CUMULATED 
(32) 


VOLUME  COSTS 
ANNUAL  CUMULATED 
(33)  (34) 


2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 


15.20 
17.80 
16.80 
13.40 
14.  10 
17.40 
11.80 
11.10 
12.20 
12.90 


1784 
2021 
1880 
1625 
1688 
1952 
1482 
1418 
1325 
1586 


1784 
3805 
5686 
7310 
8999 
10950 
12433 
13851 
15176 
16762 


388 
392 
396 
400 
404 
408 
412 
384 
361 
424 


388 
780 
1176 
1575 
1979 
2387 
2799 
3183 
3543 
3968 


150 
151 
151 
155 
156 
158 
160 
162 
140 
165 


150 
301 
452 
607 
763 
921 
1081 
1242 
1383 
1548 


2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 


8.30 
9.00 
10.90 
13.90 
13.10 
11.90 
12.70 
15.70 
13.60 


1566 
1469 
1512 
1661 
2192 
2021 
1862 
2041 
2346 
2016 


18329 
19798 
21310 
22971 
25163 
27184 
29046 
31087 
33433 
35449 


395 
433 
403 
407 
446 
416 
420 
460 
428 
433 


4363 
4796 
5199 
5607 
6053 
6468 
6888 
7348 
7776 


207 
232 
225 
214 
240 
233 
236 
248 
240 
236 


1755 
1987 
2213 
2427 
2666 
2899 
3135 
3382 
3622 
3859 


2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 


12.10 
15.20 
16.  10 
16.70 
19.60 
18.50 
14.70 
15.50 
17.10 
13.00 


1904 
2339 
2471 
2480 
2726 
2805 
2235 
2253 
2617 
2027 


37353 
39692 
42163 
44642 
47368 
50173 
52408 
54661 
57278 
59305 


473 
441 
483 
419 
455 
498 
464 
468 
513 
511 


8682 
9123 
9606 
10025 
10479 
10977 
11440 
11909 
12421 
12933 


246 
258 
261 
254 
245 
270 
262 
255 
278 
270 


4104 
4363 
4623 
4877 
5123 
5392 
5654 
5909 
6187 
6458 
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PAGE  TYPE  5  (continued)  batch    test  problem 

TEST  1 

GAME      EQUAL   AREAS   CUT  ANNUALLY 
MANAGED, THINNED  AGE  30 

TOTAL   COST  NET    INCOME  CURRENT   VALUE  TOTAL 

YEAR  ANNUAL  CUMULATED  ANNUAL  CUMULATED  GROWING   STOCK  NET  WORTH 

(35)  (36)  (37)  (38)  (39)  (40) 


1  538  538  1247                 1247  64099  65346 

2  543  1081  1477                 2724  75570  78294 

3  547  1628  1334                 4058  71139  75197 

4  555  2182  1070                 5128  55955  61083 

5  560  2742  1128                 6256  58989  65246 

6  565  3308  1386                 7643  73641  81284 

7  572  3880  910                 8553  48657  57210 

8  545  4425  873                  9426  45464  54890 

9  501  4926  824  10250  50507  60757 
1C                       590  5516  997               11247  53551  64798 

11  540  6056  659               11905  41104  53009 

12  602  6658  565               12470  32952  45422 

13  574  7232  657               13127  36007  49134 

14  557  7789  708               13835  44593  58428 

15  621  8410  1065               14900  58001  72901 

16  592  9002  1020               15920  54315  70235 

17  589  9591  840               16760  48919  65679 

18  640  10231  939               17699  52431  70130 

19  611  10841  1246               18945  65832  84777 

20  600  11442  920               19865  56440  76305 

21  660  12102  867               20732  49565  70296 

22  630  12731  1156               21887  63512  85399 

23  673  13405  1212               23099  67506  90605 

24  611  14016  1347                24446  70106  94552 

25  628  14644  1384               25830  83350  109180 

26  695  15339  1437               27267  78246  105513 

27  663  16002  1114               28380  60917  89298 

28  649  16651  1039                29420  64611  94031 

29  716  17367  1279               30698  71762  102461 

30  717  18084  904               31603  53038  84641 

PAGE  TYPE  6  PRESENT  WORTH  AND  RATE  EARNED 

BATCH      TEST  PROBLEM 
TEST  1 

GAME      EQUAL   AREAS  CUT  ANNUALLY 
MANAGED, THINNED  AGE  30 
YEARS   IN  PERIOD  30 

VALUE   OF    INITIAL   GROWING   STOCK — $  60992.01 

VALUES  DISCOUNTED  TO  PRESENT  (DOLLARS) 

COMPOUND                                 FUTURE  STOCK  NET 

RATE                                 GROWING  ALL  PLUS  ALL  PRESENT 

(PERCENT)                              STOCK  INCOMES  INCOMES  COSTS  WORTH 

1.0  39350.04  50396.60                           89746.64  15451.27  13303.36 

1.5  33931.66  46618.55                         80550.21  14330.44  5227.76 

2.0  29280.69  43222.05                         72502.74  13320.23  -1809.50 

2.5  25285.44  40163.29                           65448.73  12408.09  -7951.37 

3.0  21850.92  37403.89                         59254.81  11583.01  -13320.21 

3.5  18896.26  34910.23                           53806.49  10835.34  -18020.87 

4.0  16352.58  32652.78                         49005.36  10156.61  -22143.27 

4.5  14161.12  30605.60                           44766.72  9539.35  -25764.64 

5.0  12271.78  28745.87                           41017.65  8977.00  -28951.36 

5.5  10641.74  27053.48                           37695.22  8463.76  -31760.55 

6.0  9234.44  25510.68                           34745.12  7994.52  -34241.41 

6.5  8018.59  24101.82                         32120.41  7564.75  -36436.36 

7.0  6967.44  22813.03                           29780.47  7170.46  -38382.01 

7.5  6058.04  21632.05                           27690.10  6808.09  -40110.31 

8.0  5270.77  20548.01                           25818.77  6474.49  -41647.73 

8.5  4588.75  19551.25                         24140.00  6166.86  -43018.87 

9.0  3997.53  18633.19                         22630.72  5882.70  -44243.99 

9.5  3484.68  17786.21                           21270.89  5619.79  -45340.91 

10.0  3039.53  17003.50                         20043.02  5376.15  -46325.14 

10.5  2652.88  16278.99                         18931.87  5150.01  -47210.15 
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PAGE  TYPE  3 


ALTERNATIVES   FOR    THIS  GAME 
BATCH      TEST  PROBLEM 
TEST  1 

GAME  VARY  CUT  WITH  PRICE 
MANAGED, THINNED   AGE  30 


NUMBER   OF   YEARS   PER   GAME  30 


CRITICAL  PRICES 
ALLOWABLE  CUT 
MINIMUM   CUTTING  AGE 


12.00 
5 

130 


99.00 
10 
120 


ACRES    IN   WORKING  CIRCLE  915 

MINIMUM  VALUES  FOR    INCLUSION   IN  TOTALS 

AGE,   FOR  GROWING   STOCK  40 

M   BD.    FT.,    FOR   GROWING   STOCK  1.5 

CU.   FT.,    FOR  COMMERCIAL  CUT  400 

M  BD.   FT.,   FOR  COMMERCIAL   CUT  3.0 

M   BD.    FT. ,    FOR    SALVAGE  1.5 


COSTS   IN  FIRST  YEAR  OF  GAME 

PER   ACRE    (ANNUAL)  .20 

PER   100  CU.    FT.    HARVESTED  .05 

PER  M  BD.   FT.  1.56 

THIN  ONE   ACRE  25.00 

PLANT   ONE   ACRE  30.00 

CLEANUP   OF  ONE  ACRE  25.00 

RATE   OF    INCREASE    IN  COSTS  .01 


ACRES    PLANTED  ANNUALLY 
PERCENT  OF   ACRES   LOST  ANNUALLY 
M   BD.    FT.    IN  SHELTERWOOD 
REGENERATION  PERIOD 


RELATIVE   VALUE   OF    INTERMEDIATE  CUTS 

STUMPAGE  PRICE,  CU.  FT.  1.00 
STUMPAGE   PRICE,    BD.    FT.  .85 


PSEUDORANCOM   NUMBER  GENERATOR 


21.0 
2222.0 


PAGE  TYPE  4 


INITIAL   DISTRIBUTION  OF   ACRES  BY  AGE 
BATCH     TEST  PROBLEM 
TEST  1 

GAME     VARY  CUT  WITH  PRICE 
MANAGED, THINNED  AGE  30 
YEAR  WITHIN  GAME  0 


AGE(DECADE) 


AGE(YEAR) 
V  5 
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PAGE  TYPE  4  (continued) 


GAME      VARY   CUT   WITH  PRICE 
MANAGED, THINNED   AGE  30 
YEAR   WITHIN  GAME  30 


AGE  I  DECADE ) 


AGE  I  YEAR ) 
V  5 


PAGE  TYPE  5  BATCH     TEST  PROBLEM 

TEST  1 

GAME  VARY  CUT  WITH  PRICE 
MANAGED, THINNED   AGE  30 


ALL JWABLE 


ACTUAL  CUT 


CUMUL  CUT 
CU.FT.  MBF 
(5)  (6) 


GRSTK  VOL 
CU.FT.  > 
(7)  I 


TOTAL  VOL 


CU.FT. 
19) 


MBF 
(10) 


120 
120 
120 
130 
130 
120 
130 
130 
130 
130 


16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
14719 
16135 


129 
128 
126 
0 
78 
129 


16135 
32270 
48405 
64540 
80675 
96810 
112945 
129080 
143799 
159934 


129 
257 
383 
383 
461 
590 
655 
720 
798 
890 


145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 


4389 
4339 
4294 
4410 
4426 
4380 
4413 
4446 
4468 
4467 


161243 
177378 
193513 
209648 
225783 
241918 
258053 
274188 
288907 
305042 


4518 
4596 
4677 
4793 
4887 
4970 
5068 
5166 
5266 
5357 


130 
130 
130 
130 
130 
130 
130 
130 
120 
130 


16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
16135 
15450 


103 
116 
112 

54 
125 
140 

87 
120 
166 
121 


176069 
192204 
208339 
224474 
240609 
256744 
272879 
289014 
305149 
320599 


993 
1109 
1221 
1275 
1400 
1540 
1627 
1747 
1913 
2034 


145108 
145108 
145108 
145108 
145108 
145108 
145108 
145108 
144424 
144354 


4517 
4551 
4585 
4636 
4636 
4635 
4679 
4678 
4625 
4638 


321177 
337312 
353447 
369582 
385717 
401852 
417987 
434122 
449573 
464953 


5510 
5660 
5806 
5911 
6036 
6175 
6306 
6425 
6538 
6672 


130 
120 
120 
120 
120 
120 
130 
120 
120 
130 


16135 
16135 
15427 
16135 
16135 
16135 
16135 
16135 
15708 
16135 


115 
160 
159 
153 
156 
166 
111 
155 
178 


336734 
352869 
368296 
384431 
400566 
416701 
432836 
448971 
464679 
480814 


2149 
2309 
2468 
2621 
2777 
2943 
3054 
3209 
3387 
3430 


144284 
144215 
144145 
144075 
144006 
143936 
143866 
143797 
144154 
144094 


4637 
4588 
4537 
4484 
4446 
4393 
4392 
4353 
4302 
4401 


481018 
497084 
512441 
528506 
544572 
560637 
576702 
592768 
608833 
624908 


6786 
6897 
7005 
7105 
7223 
7336 
7446 
7562 
7689 
7831 
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PAGE  TYPE  5  (continued)  batch    test  problem 

TEST  1 

GAME  VARY  CUT  WITH  PRICE 
MANAGED, THINNED   AGE  30 


0-9 

( 12 : 


10-19 
(13) 


30-39 
115) 


50-59 
(17) 


AGE  CLASSES 
60-69      70-79  80-89 
(18)         (19)  (20) 


90-99 
(21  ) 


140-179 
(  26) 


BATCH      TEST  PROBLEM 
TEST  1 

GAME  VARY  CUT  WITH  PRICE 
MANAGED, THINNED   AGE  30 


STUMP AGE  PRICE 
100  CU.FT.  MBF 
(27)  ( 28 ) 


STUMPAGE  INCOME 
ANNUAL  CUMULATED 
(29)  (30) 


AREA  COSTS 
ANNUAL  CUMULATED 
(31)  (32) 


VOLUME  COSTS 
ANNUAL  CUMULATED 
(33)  (34) 


2.50 
2.50 
2.50 
2.50 
2.50 
2. 50 
2.50 
2.50 
2.50 
2.50 


15.20 
17.80 
16.80 
13.40 
14.10 
17.40 
11.80 
11.10 
12.20 
12.90 


2363 
2682 
2489 
403 
1501 
2614 
1169 
1124 
1325 
1586 


2363 
5044 
7533 
7936 
9438 
12052 
13221 
14345 
15670 
17256 


392 
396 
400 
404 
408 
412 
384 
361 
424 


388 
780 
1176 
1575 
1979 
2387 
2799 
3183 
3543 
3968 


209 
210 
209 

135 
220 
116 
117 
140 
165 


209 
419 
628 
636 
771 
991 
1  107 
1224 
1364 
1530 


2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2. 50 
2.50 
2.50 


10.  10 
8.  30 
9.00 
10.90 
13.90 
13.  10 
11.90 
12.  70 
15.  70 
13.60 


1462 
1385 
1424 
991 
2151 
2255 
1439 
1945 
3024 
2039 


18718 
20103 
21527 
22518 
24669 
26924 
28363 
30308 
33333 
35372 


395 
433 
403 
407 
446 
416 
420 
460 
428 
433 


4363 
4796 
5199 

5607 
6053 
6468 


189 

215 
208 
105 
235 
265 
171 
234 
321 
239 


1719 
1933 
2141 
2246 
2481 
2746 
2916 
3150 
3471 
3710 


2.50 
2.50 
2.50 
2.50 
2.5C 
2.50 
2.50 
2.50 


12.10 
15.20 
16.10 
16.70 
19.60 
18.50 
14.70 
15.50 
17.10 
13.00 


1802 
2815 
2962 
2970 
3484 
3499 
2049 
2821 
3470 
978 


37174 
39989 
42952 
45922 
49406 
52905 
54954 
57774 
61245 
62223 


473 
441 
483 
419 
455 
498 
464 
468 
513 
511 


9123 
9606 
10025 
10479 
10977 
11440 
11909 
12421 
12933 


230 
319 
320 
311 
322 
345 
236 
330 
381 
102 


3940 
4259 
4579 
4890 
5212 
5557 
5794 
6123 
6504 
6606 
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PAGE  TYPE  5  (continued) 


LAT'H      TFSI  PROBLEM 
TEST  1 
GAME  VARY 
MANAGED, THI 


TOTAL  COST 
ANNUAL  CUMULATE!"! 
(35)  136) 


MFT  INCOME 
>iN'JAL  CUHUL  ATEU 

I3M  138) 


CURRENT  VALUE 
CROWING  "STOCK 
(39) 


TOTAL 
NET  WORTH 
CO) 


597 
602 
604 
408 
539 
628 
528 
501 
501 
590 


597 
1  199 
1803 
2211 
2750 
3378 
3905 
4406 
4908 
5497 


1765 

2080 
1885 
-5 
963 
1986 
641 
623 
824 
997 


1765 
3845 
5730 
5725 
6688 
8674 
9316 
9938 
10762 
11759 


63417 
73946 
68856 
55558 
58777 
72580 
48313 
45464 
50507 
53551 


65182 
77792 
74586 
61283 
65465 
81255 
57628 
55403 
61269 
65310 


495 
557 
530 
512 
625 
596 
543 
645 
691 
603 


5993 
6550 
7079 
7592 
8216 
8813 
9356 
10001 
10692 
11296 


1093 
1046 
587 
970 
1843 
940 


12200 
12597 
13070 
13548 
14641 
15687 
16273 
17243 
19086 
20026 


41389 
33399 
36751 
45848 
59666 
55867 
50660 
54307 
67500 
57855 


53589 
45996 
49821 
59396 
74307 
71554 
66934 
71550 
86586 
77882 


664 
710 
753 
689 
705 
770 
658 
724 
787 
549 


1  1960 
12670 
13423 
14112 
14817 
15587 
16245 
16969 
17756 
18305 


887 
1710 
1791 
1934 
2065 
2056 
1085 
1533 
1  794 
23 


20913 
22624 
24415 
26348 
28413 
30469 
31554 
33087 
34861 
34904 


50797 
64429 
67750 
69601 
81844 
75995 
59188 
62090 
68202 
51552 


71710 
87052 
92165 
95949 
1 10258 
106464 
90742 
95177 
103083 
86456 


PAGE  TYPE  6 


PRESENT   WORTH   AND  RATE  EARNED 
BATCH      TEST  PROBLEM 
TEST  1 

GAME      VARY  CUT  WITH  PRICE 
MANAGED, THINNED  AGE  30 
YEARS    IN   PERIOD  30 


VALUE   OF    INITIAL   GROWING  STOCK- 


VALUES   DISCOUNTED   TO   PRESENT  (DOLLARS) 


COMPOUND 

RATE 
(PERCENT) 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

6.0 

6.5 

7.0 

7.5 

8.0 

8.5 

9.0 

9.5 
10.0 
10.5 


FUTURE 
GROWING 
STOCK 

38247.39 

32980.84 

28460. 19 

24576.89 

21238.62 

18366. 76 

15894.35 

13764.30 

11927.90 

10343.54 

8975.67 

7793.90 

6772.20 

5888.29 

5123.07 

4460. 16 

3885.51 

3387.03 

2954.35 

2578.54 


ALL 
INCOMES 

52747.65 

48739.60 

45142.72 

41909.46 

38998. 17 

36372.40 

34000.11 

31853.20 

29906.94 

28139.54 

26531.81 

25066.81 

23729. 57 

22506.84 

21386.90 

20359.34 

19414.93 

18545.45 

17743.61 

17002.89 


STOCK 
PLUS 
INCOMES 

90995.04 

81720.44 

73602.91 

66486. 35 

60236.79 

54739. 15 

49894.46 

45617.51 

41834.84 

38483.08 

35507.48 

32860.71 

30501.77 

28395. 13 

26509.97 

24819.50 

23300.44 

21932.48 

20697.96 

19581.44 


ALL 
COSTS 

15614.36 

14470.52 

13440.43 

12511. 14 

11671.31 

10910.99 

10221.45 

9595.00 

9024.86 

8505.07 

8030.34 

7596. 04 

7198.01 

6832.62 

6496. 61 

6187.10 

5901. 53 

5637.60 

5393.28 

5166. 75 


NET 
PRESENT 
WORTH 

14388. 67 

6257.91 

-829.53 

-7016.80 

-12426.53 

-17163.85 

-21319.00 

-24969.50 

-28182.04 

-31014.00 

-33514.87 

-35727. 34 

-37688.26 

-39429.50 

-40978.65 

-42359.61 

-43593. 10 

-44697. 13 

-45687.33 

-46577.33 


-  31  - 


PAGE   TYPE   7  COMPARISON  OF  ALTERNATIVES 

BATCH     TEST  PROBLEM 
TEST  1 

MANAGED, THINNED   AGE  30 
COLUMN  10 


4530 
4621 
4716 
4807 
4897 


4518 
4596 
4677 
4793 
4887 


4992 
5083 
5173 
5273 
5364 


4970 
5068 
5166 
5266 
5357 


5364 
6632 
7893 


5357 
6672 
7831 
0 


90 
100 


110 
120 
130 
140 
150 


COMPARISON  OF  ALTERNATIVES 
BATCH     TEST  PROBLEM 
TEST  1 

MANAGED, THINNED   AGE  30 
COLUMN  40 


65346 
78294 
75197 
61083 
65246 


65182 
77792 
74586 
61283 
65465 


81284 
57210 
54890 
60757 
64798 


81255 
57628 
55403 
61269 
65310 


64798 
76305 
84641 


65310 
77882 
86456 


90 
100 


110 
120 
130 
140 
150 
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